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We study the evolution of phase-transition-generated cosmic magnetic fields coupled to the 
primeval cosmic plasma in turbulent and viscous free-streaming regimes. The evolution laws for 
the magnetic energy density and correlation length, both in helical and non-helical cases, are found 
by solving the autoinduction and Navier-Stokes equations in mean-field approximation. Analyt- 
ical results are derived in Minkowski spacetime and then extended to the case of a Friedmann 
universe with zero spatial curvature, both in radiation and matter dominated eras. The three pos- 
sible viscous free-streaming phases are characterized by a drag term in the Navier-Stokes equation 
j.-^ I which depends on the free-steaming properties of neutrinos, photons, or hydrogen atoms, respec- 

»vj . tively. In the case of non-helical magnetic fields, the magnetic intensity B and the magnetic corre- 

lation length ^B evolve asymptotically with the temperature T as B{T) ~ usiNiViY'^ [T /TtY^ and 
(,b{T) ~ K^{NiVi)^^ {T /Ti)^'^ . Here, Ti, Ni and Vi are, respectively, the temperature, the number 
^^ , of magnetic domains per horizon, and the bulk velocity at the onset of the particular regime. The 

'^N ' coefficients kb, kj, qi, Q2, Qs, and g4, depend on the index of the assumed initial power-law magnetic 

l^-s ' spectrum, p, and on the particular regime, with the order-one constants kb and k^ depending also 

on the cut-off adopted for the initial magnetic spectrum. In the helical case, the quasi-conservation 
of the magnetic helicity implies, apart from logarithmic corrections and a factor proportional to the 
initial fractional helicity, power-like evolution laws equal to those in the non-helical case, but with 
p equal to zero. 
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I. INTRODUCTION 



o 
u 

o 

^ [ All observed galaxies and cluster of galaxies present microgauss, large-scale magnetic fields. The origin of these 

I I - cosmic magnetic fields has not yet been fully understood, although several proposals, especially in the last years, have 

[ been put forward to explain why our universe is magnetized (for reviews on cosmic magnetic fields, see [l|-l5|). These 

proposals belong to two distinct classes of generating mechanisms characterized, roughly speaking, by the time when 

they operate. Astrophysical mechanisms work on at the epoch of large-scale structure formation or later, and can 

be based on the Biermann battery effect [6[ in the first supernova remnants [7|, or on the Weibel instability [8| in 

f^ ' intergalactic plasmas [9|. Mechanisms performing in the early universe can, on they part, classified in mechanisms 

^^ I operating in the very early universe, namely during an infiationary epoch of the universe [lO[, and mechanisms 



tH- [ operating after inflation |ll| - |l4| . for example during electroweak (QED) or quark- hadron (QCD) cosmological phase 
^^ . transitions |15j . 

C^ ■ Magnetic fields generated during large-scale structure formation through astrophysical processes usually have large 

\ I . correlation scales but low intensities. Then, the precondition to explain the observed fields is an amplification of 
^ ' these seeds. In principle, dynamo mechanisms operating in gravitationally bound large-scale structures (such as 
J^ I galaxies and clusters of galaxies) could substantially increase the intensity of such fields up to observed values |16l |. 
\^ ■ However, astrophysical mechanisms cannot explain the presence of large-scale magnetic fields recently detected in 
Cd , cosmic voids |17| . 

Inflation-produced magnetic flelds have the potentiality to explain the origin of cosmic magnetic fields since they 
possess a large correlation. The unpleasant feature of this class of mechanisms is that, in order to obtain astrophysically 
interesting magnetic intensities, one has to introduce some nonstandard (to wit, speculative) interaction term in 



the photon field Lagrangian. A possible exception is represented by a generating mechanism first noticed in 18 1 
and successively analyzed in the series of papers [19|. Here it is shown that, starting from the standard Maxwell 
Lagrangian and assuming a small negative curvature term (not present in all the other infiationary mechanisms), a 
very strong magnetic field can emerge as the result of photons creation from a varying gravitational field. The validity 
of the results in theses works, however, has been put into question in [20| and [2l|. According to the authors of [21| 
indeed, the intensity of the produced field would be today, in the best case scenario, as low as Bq ~ 10~^^G and then 
astrophysically unimportant. Another exception (the last one to our knowledge) to the use of nonstandard physics 



for the case of inflationary magnetic fields is represented by the recent work [22|. Here, criticisms to all previous 
inflationary generating mechanisms for cosmic magnetic fields have been drawn. The main criticism concerns the use 
of unrenormalized, and then unphysical, vacuum magnetic fluctuations instead of the renormalized ones (the only 



ones which are physically acceptable). In [22|, it is shown that in standard Maxwell theory, the actual magnetic 
field resulting from inflationary renormalized vacuum magnetic fluctuations is scale independent and has an intensity 
which depends only on the scale of inflation. If this scale is around lO-'^^GeV (a plausible value for such a scale), that 
field directly account for galactic and galaxy cluster magnetic fields. 

Magnetic fields created during cosmological phase transitions have small length scales. This is due to the fact that 
microphysical processes which participate to the generation of the magnetic field are necessarily uncorrelated on scale 
greater than the Hubble scale which represent an (event) horizon for all causal physical processes. Nevertheless, the 
primordial universe could had been very turbulent at the time of QED and QCD phase transitions. This means 
that any magnetic field coupled to the primeval plasma could had been processed by magnetohydrodynamic (MHD) 
effects. Indeed, it is well-known in the literature that turbulence effects increase the typical magnetic length scale ^]. 
Moreover, if the magnetic field is helical the growth of magnetic correlation could be much more significant. It is 
interesting to observe that if the large-scale magnetic fields we observe today were relic helical fields from cosmological 
phase transitions (or inflation), then this would be a manifestation of a macroscopic and primeval P and CP violation, 
since magnetic helicity is odd under discrete P and CP transformations. ^ Moreover, models for generating helical 
magnetic fields in the early universe there exist in the literature 2J, |25| . 



It is worth noting that magnetic fields produced in the early universe are subject to a variety of constraints coming 
from the fact that the presence of a primordial magnetic field could spoil the predictions of the standard cosmological 



model. The most important limits came from the study of Big Bang Nucleosynthesis [26| and Cosmic Microwave 



Background (CMB) anisotropics [27l- l29| . In particular, large-scale, primordial helical magnetic fields could leave 
peculiar signatures in CMB radiation [30|, since they would introduce a parity-odd cross-correlations of the CMB 
anisotropics, not present in the case of a non-helical magnetic field. 

Also, primordial helical magnetic fields could affect the phenomenology of axions since the latter are coupled, 
through the axion-photon interaction term, to external magnetic fields. Indeed, the generation of helical magnetic 
fields before the axion coherent oscillations start at a temperature around few GeV, could be in contradiction with 
the existence of the axion [3l| , since a strong helicity production would in turn produce too nuich of an axion relic 
abundance, in disagreement with astrophysical observations. 

Finally, it is interesting to observe that the presence of large-scale helical magnetic fields in the actual universe 



could be directly detected by analyzing the propagation properties of charged cosmic rays, as shown in 32 1. 

Magnetic fields generated during inflation or cosmological phase transitions can evolve as the universe expands, since 
their properties can be modified by magnetohydrodynamic turbulent effects due to the coupling with the primeval 
cosmic plasma. As recently shown in [33| . typical inflation-produced magnetic fields remain almost unchanged on 
scales of cosmological interest even in the presence of a turbulent plasma (as in the case, for example, of a cosmic 
phase transition) although on small scales, which are however astrophysically uninteresting, their power get gradually 
suppressed. 

The fact that the universe during a cosmological phase transition could be very turbulent was first realized in j34| , 
where the importance of MHD turbulent effects on the evolution of phase-transition-produced magnetic fields was 
stressed. Here, in order to study the correlation properties of an evolving cosmic magnetic field, however, a simple 
toy model was employed, which replaced the full MHD equations. Since then, there are been many efforts to study. 



more accurately, the evolution of turbulent magnetic fields |3 51451 1 . 

However, only in the work 52], it was realized that, other than turbulent phases, phase-transition-generated mag- 
netic fields can experience dissipation processes induced by free-streaming neutrinos, photons, or hydrogen atoms. 
Indeed, such cosmic components can have mean free paths which exceed the typical correlation scale of a primordial 
magnetic field, so can be considered as free-streaming species with respect to the magnetic field itself. In this case, the 
dissipation of magnetic energy is ruled by a "drag" term in the Navier-Stokes equation, instead of the usual diffusion 
term proportional to the kinematic viscosity. 

Such a possible phase in MHD, called "dragged magnetohydrodynamics" 53J,|5J1, has been fully analyzed in J53l,l54{- 
where a systematic study of the evolution of cosmic magnetic fields has been performed by direct numerical integration 
of MHD equations. To our knowledge, numerical simulations of dragged MHD have so far been performed only 
in ISSj, [54| where, moreover, a justification of the results has been given by using simple scaling arguments. 

The aim of this paper is to analyze the evolution of phase-transition-generated magnetic fields, in turbulent and 
viscous free-streaming regimes, by using suitable mean-field-approximated magnetohydrodynamic equations. 



^ Helical magnetic fields are indeed characterized by the fact of possessing an asymmetry between the number of left-handed and right- 
handed photon helicity states (see, e.g., \24\). 



In particular, we will derive the evolution laws for the magnetic energy and correlation length both in non-helical 



and helical cases. We anticipate that our results agree with the numerical results in J53l.l54|. 

The paper is organized as follows. In the next section, we introduce the relevant equations and the physical quantities 
of dragged MHD. Also, we find some new exact results in dragged MHD, and we introduce the transformations 
connecting the evolution of a primordial magnetic field in a (static) Minkowski spacetime to that in a flat Friedmann 
(expanding) universe. In section III, we apply some rigorous scaling arguments, first introduced by Olesen for the 
study of turbulent MHD, to dragged MHD. In section IV, we find the equations governing the evolution of the 
magnetic energy and helicity by using three different mean-field-theory approaches. In section V, we solve these 
equations for some specific initial conditions, and find the evolution laws of the magnetic intensity and correlation 
length in radiation dominated universe for the case of free-streaming neutrinos and photons. In section VI, the case 
of free-streaming photons and hydrogen atoms in matter dominated universe is analyzed. In section VII, we justify 



some scaling arguments of [53|, [5^ in the light of our exact analytical results obtained in section II. In section VIII, 
we discuss the evolution of freely-decaying turbulent magnetic fields, comparing the various and often different results 
present in the literature. Also, we complete our work started in j48.] by extending our previous results, obtained in 
Minkowski spacetime, to the case of a fiat Friedmann universe. In section IX, we study the dependence of our results 
on the choice of the initial cut-off of the assumed, initial magnetic power spectrum. In section X, we discuss our results 



and compare them to the results of [53j, |5J|. In section XI, we draw our conclusions and give possible suggestions for 



further investigations. Finally, appendixes A,B, and C contain some technical details of our analytical computations. 

II. DRAGGED MAGNETOHYDRODYNAMICS 

In this section, we introduce the equations and physical quantities of interest for dragged magnetohydrodynamics. 
Some new exact results are derived, some scaling arguments are discussed, and the transformations relating the 
evolution of a magnetic field in Minkowski and Friedmann spacetimes are introduced. 

Ila. Preliminaries 

We start by writing down the incompressible ^ magnetohydrodynamic equations in the presence of a drag term and 
in Minkowski spacetime 52|, 



Dv 

-— == (V X va) X va - av, (1) 

at 



dt 



V X (v X v^) + ?7V'v^, (2) 



where D/di = dt + v-V is the so-called hydrodynamical derivative, v is the velocity of bulk fiuid motion and V • v = 0, 
VA = B/^/p -I- P is the Alfven velocity, B is the magnetic field, and p and P are the energy density and the thermal 
pressure of the fiuid, respectively. The frictional coefficient a and the resistivity rj are dissipative parameters, and are 
determined by microscopic physics. 

We will consider throughout the paper the case of constant p and P. In the case of a Friedmann universe, where p 
and P evolve in time, the quantities of interest, as we will see in section lib, are p and P, instead of p and P, with 
the "tilded" quantities being time independent. For this reason, and for the sake of simplicity, we can take p + P = 1 
in Eq. ([1]), yielding v^i = B. 

The Navier-Stokes equation H]) and the autoinduction equation ([2]) can then be written, respectively, as 

Dv 

^^F.-F„ (3) 

— =Vx(vxB)+r?V2B, (4) 



^ The kinetic energy associated to a turbulent fluid is E^ = {l/2)pv'^, where p is the energy density of the universe, and v is the typical 
velocity of bulk fluid motion associated to turbulence. Incompressibility of the primordial plasma means that the acoustic Mach number, 
Ms = v/vb, where Vb is the speed of sound, is much less then unity [23| . For example, in the radiation era, when the cosmic fluid is 
ultra-relativistic and then Vb = I/a/S, the incompressibility condition reads v <C l/VS- 



where we have introduced the Lorentz force F^ = J x B, and the damping force F^ = av, with J = V x B being 
the magnetic current. Let us now introduce the relevant physical quantities in (dragged) MHD. The magnetic and 
kinetic energy densities are defined as 

Esit) = i /"d^xB2(x,i) = / dk£Bik,t), (5) 

Ev{t) = i fd'^xv\^,t) = / dk£,{k,t), (6) 

where £b — 27rfc^ |B(k)p and £„ — 27rfc^ |v(k)p are the magnetic and kinetic energy density spectra. Here, B(k) and 
v(k) are the magnetic and kinetic fields in Fourier space, with k being the wavenumber and k = |k|. The magnetic 
and cross helicity densities are 

//•oo 
d^xA-B= / dknB{k,t), (7) 

Hc{t)^ jd^xv-B= I dknc{k,t), (8) 

respectively, where Hb = 47rfc^A(k)-B*(k) and He — 47r/c^v(k) •B*(k) are the magnetic and cross helicity density 
spectra, A is the vector potential, and an asterisk denotes complex conjugation. It is worth noticing that, for all 
magnetic and kinetic field configurations, the magnetic and cross helicity spectra satisfy the following realizability 
conditions: IV-b] < ik'^^B and |Hc| < 2(£„fs)^^^- The total energy Etot = Eb + Ey and the cross-helicity are 
conserved quantities when a = 77 = 0, while the magnetic helicity is conserved when rj — 0. This follows directly from 
their evolution laws, which can be straightforwardly derived from the MHD equations ^ and (jlj. They are 

^^-2aEy-r,jd\3\ (9) 

-^ = -aHe-r, d'^xJ-u,, (10) 

^ = ^2rjJd'x3-B, (11) 

where a; = V x v is the so-called vorticity. We also introduce the relevant length scales in (dragged) MHD, i.e., 
the so-called correlation lengths, which are the characteristic lengths associated with the large magnetic and kinetic 
energy eddies of turbulence. These are defined by 

r (t) lo^dkk-^^BJk) 

^^^'^- J.-'dkSBik) ' ^''^ 

_ J-dkk-^£Ak) 
^^^'>- !^dk£y{k) • ^'^^ 

With the aid of the magnetic and kinetic correlation lengths, it is straightforward to transform the above "local" 
realizability conditions on the magnetic and cross helicity spectra, into the following "global" realizability conditions: 
\E[b\ < '^S.bEb and \Hc\ < min[£'tot, 2 jdk [EyEBY^"^ ]• Finally, it is useful to define the kinetic Reynolds number, the 
magnetic Reynolds number, and the Prandtl number as Re = v/{la), Re^ = vl/rj^ and Pr — Re^/Re, respectively, 
where v and I are the typical velocity and length scale of the fluid motion. 

lib. The equilibrium state 



A dragged phase develops when the kinetic Reynolds number is small [531 . I54| . For this reason, we work in the case 
where Re ^ 1. Defining the "drag time" r^ = a"^, we have |Dtv|/|Frf| ~ Td/reddy = Re, where we used dt ~ l/roddy 
and introduced the so-called (kinetic) eddy turnover time Toddy — l/v- Accordingly, for very low kinetic Reynolds 
numbers, we can neglect the left-hand side of Navier-Stokes equation ([3|). This means that the system will approach 
asymptotically a state where the Lorentz force equilibrates the damping force 53, 54], 

FL=Frf. (14) 



In this "equilibrium state" , the autoinduction equation Q reduces to an equation which depends only on the magnetic 
field, 



dB 

'dt 



= rrfV X [B X (B X V X B)] + t^V^B, 



(15) 



and will be basis for the study of the evolution of the magnetic field itself. 

Multiplying both sides of Eq. ([T^ by B, wc find that the cross helicity is zero in the equilibrium state, 



H,{t) = 0. (16) 

Also, there is no equipartition between kinetic and magnetic energy in this state. Indeed, defining the F-ratio 



we get 



Eb 



r - Re < 1, 



(17) 



(18) 



since V ~ v^ /B^ ^ Re |Fi|/|F(j| = Re, where B is the typical intensity of the magnetic field. Moreover, Td and Toddy 
are related by 



Td ^ rToddyi 



(19) 



since 1~ |Fd|/|rL| - Freddy /t^ . Let us observe that \E^\/\Eb\ -- Re^/^ \aE^\/\EB\ ^Re^'"^ , &nd\T]jd'^x3'^\/\EB\ - 
Re ' Re^ , where a dot indicates a time derivative. Consequently, in the limit of low kinetic Reynolds number and 
large magnetic Reynolds number, Eq. Q can be approximated with 



dEB 
dt 



-2aEy. 



(20) 



The above equation connects the evolutions of magnetic and kinetic energies, and will be useful in the following. 
Finally, since we are considering a very large Prandtl number, we can neglect the second term in the right-hand side 
of Eq. (fTO| . In fact, we have \ri Jd^xJ ■ (jLj\/\aHc\ ^ 1/Pr. Accordingly, we get He — —ctHc, whose solution 



H,{t) = H,{U)cicp 



- / dt'a{t') 

ti 



(21) 



gives the evolution law of the cross helicity in the equilibrium state. In the next subsection, we show that the above 
equation is indeed (asymptotically) compatible with Eq. (J16p . 



lie. Scaling arguments 



Let us suppose, and this is indeed the case in a cosmological context 53l|54|, that the drag coefficient a evolves in 
time as a simple power law. 



(22) 



a{t) = a{ti) ( — 



where ti is a reference time. 

In numerical simulation of dragged magnetohydrodynamics 53l IsJ] , it is observed that the magnetic energy scales 
in time as a power law if the initial magnetic spectrum is assumed to be a simple power of the wavenumber. Let us 
then parameterize the magnetic energy as 



EB{t)=EB{t,) 



(23) 



Inserting Eqs. ((22|) and (|23p into Eq. (l20l) . we find the evolution laws for the kinetic energy and F-ratio, 



^x-(l+a+/3) 



r{t)^T{u)[i 



-(1+q) 



(24) 
(25) 



respectively, where T{ti) = {(3/2)[Td{ti)/ti]. From this last relation we get /3 > 0, which means that the magnetic 
energy decreases in time during a dragged MHD phase (as in the case of turbulent MHD). 

1 /2 

Moreover, remembering that F ^ Re — v/al and taking v ^ Ev , we find that the typical length scale evolves in 
time as 



giving, in turn, the evolution law 



lit) (X i(l--«/2, 



Teddy (i) « t 



for the eddy turnover time. Finally, inserting Eq. ([2^ in Eq. (I^TI) . we get 



H,{t)^H,{t,)exp\ 



1 



if a 7^ —1, and 



H,{t) =. HS 



TdiU) l + a 



l+a 



(26) 

(27) 

(28) 
(29) 



if a = —1. In both cases, since ti ~ Teddy(ii) 2> Td{ti) [see Eqs. p^ and ([T^ ]. we get that the cross helicity goes to 
zero asymptotically {t '^ ti), in agreement with Eq. ((T6)) . 

It is amazing to observe that, even without knowing the value of the exponent /3, simple scaling arguments, when 
combined with the exact laws (|20|) and (f2T|) . give the evolution laws for the F-ratio, the eddy turnover time, and the 
cross helicity. 



lid. Expanding universe 



In the case of an expanding universe (with zero spatial curvature) described by a Friedmann- Robertson- Walker line 
element, it has been shown that the dragged MHD equations are the same as the Eqs. ([T|) and ([2]) provided that time, 
coordinates and dynamical variables are replaced by suitable quantities. The form of these new "tilded" quantities is 
different in radiation an matter eras. 



Radiation era: comoving variables- In a radiation-dominated universe the tilded quantities are 53|, l54| : 



t- 

B 

P- 
V' 



i = R ^dt, X — ;► X = Rx, 



> B = R^B, 

P 



R^P, 
R-'v, 



V — > V = V 

p ^ p ^ 

a — > (5 = Ra, 



R^P, 



where R{t) is the expansion parameter, t is the so-called conformal time, and we note that v is not scaled. 
Matter era: supercomoving variables- In a matter-dominated universe the tilded quantities are [53|, |5J] : 



t^£=/«-3/^.<, X^i^ 



i?X, 



B -^ B = R^B, 


V — > V = 


= i?i/2v. 


P^P^ R^P, 


P ^ P - 


= R^P, 


T] ^ fi = R'^^^r], 


a — > a = 


-- R^/'a. 



(30) 

(31) 

(32) 
(33) 



(34) 

(35) 
(36) 

(37) 



Due to the formal coincidence of the dragged MHD equations in Minkowski and Friedmann spacetimes, we can analyze 
the evolution properties of a primordial magnetic field in both cases in a similar way. For definiteness, we consider first 
the case of a non-expanding universe and then translate the results to the case of interest of an expanding universe. 
We observe that, since p oa P oa R~'^ in the radiation-dominated era and p ex. P (x R~^ in the matter-dominated 
era, the quantity p -I- P is a constant during the evolution of the universe. This result was used in section Ila. 

III. OLESEN'S APPROACH 



In this section, using a scaling approach first introduced by Olesen in [35| for the case of freely-decaying MHD 
turbulence, we express the unknown exponent (3 introduced in section lie as a function of the exponents of the 
assumed initial power-law magnetic spectrum, £B{k,t = 0) oc /c^, and of the power law for the drag coefficient, 
a{t) ex r. 

Non-helical case- We are interested to the case of large magnetic Reynolds numbers, and then we neglect the 
dissipative term in the autoinduction equation Q. We can then write the dragged MHD equations as 

V = - (V X B) X B, (38) 

a 

— =Vx(vxB). (39) 

A direct inspection show that, under the scaling transformations x — )■ £x and t — )■ £^~"'t, equations (|38p and ([M)) 
admit solutions of the type 

v(^x,^i-"t)=^"v(x,t), (40) 

B(^x,^i-"i)=^(i+"+'")/2B(x,t), (41) 

a(£i-"i)=r"a(i), (42) 

where ^ > is the "scaling factor", and u and m real parameters. Differentiating Eq. P^ with respect to £, and 
putting £ = 1 afterwards, we get a{t) oc i™/(i^"), and then m — a(l — u), since we are assuming that a oc i°. Defining 
p = —{2 + u + m) and using Eqs. PO]) and (|4ip . we straightforwardly obtain 

Seik, t) = AbF il^B{k^^'^ t^^"), (43) 

f„(fc,i) = A.fc>4fc3+Pii-'^), (44) 

where we have defined r = [3(1 + a) + 2p]/{l — a) for notational convenience. Here, A^ and A^, are constants, while 
ipB and '01, are arbitrary scaling-invariant functions of their arguments. 

What is observed in dragged MHD [53|, |5J| (as well as in turbulent MHD 53[ |54|) is that the evolution of the 



magnetic spectrum in the non- helical case proceeds through the so-called "selective decay" (see section V), so that 
the initial magnetic spectrum retains its form (for lengths below the characteristic dissipation scale) for all times. 
This means that if we assume that the initial magnetic spectrum is a simple power of the wavenumber fc, then the 
coefficient p is the index of that power. In fact, we have Ss^k, 0) ~ AsfcP iPb{0), a result which is true also in turbulent 
MHD (see fsl). Inserting Eq. ^ into Eq. © and Eq. (jHl) into Eq. §^, we obtain Eqs. ^ and ^, respectively, 
where now (3 — {1 — a)(l +p)/(3 + p) is expressed as a function of p and a, 

^BCxt-'5-t-(i-'')(i+P)/(3+P), (45) 

Equations (|25|) and ([27)1 are accordingly confirmed. Also, we find that 

es ex C. oc id— «/2 = t(i-")/(3+P) ^ (47) 



confirming the scaling relation (|26|) . 

Helical case- In the case when the magnetic helicity is different from zero, we can still apply the Olesen's approach 
to find the scaling behavior of the magnetic helicity spectrum. We find, straightforwardly, 

nB{k,t)^fiBkP-'cl)B{k''+n'-''), (48) 



8 



where hb is a constant and 4>b an arbitrary scaling-invariant function. In dragged MHD [53|,|54| (as well as in turbulent 
MHD [23|, |53|, l5J|), the magnetic spectrum evolves through a phase of so-called "inverse cascade" (see section V), 
during which magnetic modes with small wavelengths get modified by the longer ones and viceversa. During this 
process, the information enclosed in the initial magnetic spectrum are lost, so that the coefficient p in Eq. (H5)) does 
not represent, necessarily, the index of the initial power-law spectrum. Indeed, since magnetic helicity turns to be 
(quasi-) conserved in dragged MHD [53[ |5J| (as well as in turbulent MHD [23|, |53|, [SJ]), we must have p = 0. In fact, 
equations ([7]) and p8)) imply Hb^I) = HB{ti){t/ti)~'-^~°-'>P^^^^P\ Accordingly, the evolution laws for the magnetic and 
kinetic energies and correlation lengths are obtained in the helical case taking p — in the evolution laws previously 
found for the non-helical case: 



Eb oc i-(i-)/3, 

E, oc i-2(2+'^)/3^ 

Equations (|45]), (gSl), (gT]), and equations (09]), ([501), dSI]), are in agreement with the results of [ssllsil. 



(49) 
(50) 
(51) 



IV. MEAN-FIELD APPROXIMATION 



In this section, working in mean-field-theory approximation, we derive the evolution integro-differential equations 
for the magnetic energy and helicity spectra, which will be solved later on in section V. We start from the "simplest" 
mean- field approximation first introduced in |55|, we then study the dragged MHD equations in the so-called "one- 
point-closure" 23 1 and "quasi- normal" [48| approximations. 



IVa. Cornwall's approximation 



It is useful to re- write the autoinduction equation (|15p as 

dB 



— = TrfV X [(J • B)B - B"J] + ryV'B. 



(52) 

;j.B) 

and B"^ ^- i (B"^), where the brackets (...) indicate a suitable average to be defined later. This operation allow us to 



We now proceed as in 55| by replacing the quadratic terms J • B and B in the above equation with J • B — > ^ 

s (...) indicate a suitable average to be defined later. This operation all 
linearize the autoinduction equation. 



dB 



Td 



— = ^ V X [(J • B)B - (B^) J] + TyV^B. 



(53) 



Whatever is the averaging operation, the averaged quantities in the above equation must satisfy a consistency relation 
coming from conservation of magnetic helicity in the case of null dissipation. Indeed, observing that Eq. (|53p gives 



OHe 



^ fd^x [(J • B)B^ - (B^)J • B] - 2ry fd^xJ ■ B, 



we get 



d-^x(3 ■B)B-' 



d^x(B'^)J-B. 



(54) 



(55) 



Assuming now, as in [55[, that the operation of averaging is just a volume average, (J • B) = Jd^xJ ■ B and 
(B^) = fd^xB^, we find that Eq. ([55)) is satisfied for all magnetic field configurations and, taking into account 
Eqs. dm and (P, that (J • B) = -HbH^vi) and (B^) = 2Eb- Finally, we can re-write Eq. ^ as 



(56) 



dB , 

— -==aBVxB + 77effV2B, 
ot 

where we have introduced the quantities asft) = —Hb Td/{6ri), ricg{t) = rj + rjT, and r]T{t) — 2EBTd/3. Equation ([56]) 
describes the well-known a— dynamo effect 23|, in which the dynamo, the turbulent diffusion, and the total effective 
diffusion coefficients, are given by aB, ?7t, and ?7cff, respectively. 



IVb. One-point-closure approximation 

Cornwall's arguments find a full justification in the framework of one-point-closure approximation, originally intro- 
duced in turbulent MHD [23|, and that we apply now to dragged MHD. 

Equation- Let us assume that v and B can be decomposed into an average part varying only on large scales and 
a weak, small-scale fluctuating part, 

v = vo + v, B = Bo + B, (57) 

with (v) = (B) = 0, and |v| ^ |vo|, |B| ^ |Bo|. In a moving coordinate system such that Vq = 0, the equation 
governing the evolution of the mean magnetic field is a dynamo equation 

-^ ^ a'sV X Bo + {vt + vW^Bo, (58) 

with dynamo and turbulent diffusion coefficients given by [23[ 

«'bW =-^ (v- V X v) ^ -^H,, (59) 

V'Tit) = ^ (V^) ^ ^ E^^^Eb Freddy. (60) 

In Eq. ([5^ . we have introduced the the so-called kinetic helicity density [23|, Hy{t) = ^d^xw ■ uj^ with u; = V x v, 
while in the last equality of Eq. ([60t . we used Eq. (flT)) . 

Starting from Eq. (|58l) and imposing the conservation of the mean magnetic helicity Jd^xAo ■ Bq — Jd^x A ■ B 
for vanishing resistivity, we get a^/d'^xBg — ri'^^Jd'^xJo ■ Bq. Since Jd'^xBg ~ Jd'^x'B — 2Eb and Jd'^x Jq • Bq — 
fd^xJ ■ B = -Hb/{2iii) [see Eq. ((TTl)], we have 2EBa'g = -[H b /{2?])] tj'j, and in turn, using Eqs. dSH) and (pOI) . 
we find Hy = {T/2'q) Hb- Therefore, we can re-write the dynamo coefficient as a'g{t) ~ — ij^ rTeddy/(6?7)- Both 
the dynamo and turbulent diffusion coefficients are then consistent with the expressions for a^ and 77T because of 
Eq. IIID. 

Solution- In order to solve the dynamo equation ([55| [or, which is the same, Eq. ([55]) ] it is useful to introduce 
the orthonormal helicity base {e_|_,e_,e3}, with e± — (ei ± ie2)/y/2 and 63 = k/fc, where {61,62,63} is a right- 
handed orthonormal base. In the helicity base, the magnetic field in Fourier space can be decomposed as B(fc,i) = 
i3+(fc,i)e+ + B-{k,t)e-, where B+{k,t) and B^{k,t) represent the positive and negative helicity components of 
B(fc,i), respectively. Equation ((56|) becomes B± = ±aBkB± — r]csk^B± in Fourier space. The solution of the above 
equation is easily found, 

B±(fc,0 = B±(fc,0)exp(±fc^,-fc24^J, (61) 

where we have defined the "dynamo" and "dissipation" lengths, ia{t) = /q dtaB and ^diss(^) — /o d.trjcs, respectively. 
In the helicity base, the energy and helicity spectra read 

fs(fc,t)=(^) (IS+P + IS-I^), (62) 

nBik,t) = ^{\B+\'-\B^\'). (63) 

To proceed further, let us suppose that \B_{k,0)\'^ = {1 - h)\B+{k,0)\'^, and define ?^g^''(fc,i) =2fc-ifB(fc,i) so that 
we can write the initial helicity spectrum as ■Hb(A:, 0) — hB^g'^^ik, 0), where Hb = h/{2 — h). Indeed, we have made 
the simplifying assumption that the initial helicity is just a fraction of the initial maximal helicity. In other words, we 
restrict our analysis only to so-called magnetic fields with (initial) "fractional helicity" . Inserting Eq. (pT|) in Eqs. (p^ 
and (p5|) we get, respectively, 

£B{k,t) = £B{k,Q) [cosh(2H„) + /iBsinh(2H„)]exp(-2fc24^J, (64) 

•HB(fc,t) =Hg^''(fc,0)[sinh(2M„) + /iBCOsh(2fc4)]exp(-2/c2£2.^J. (65) 

The above equations are integro-differential equations for the magnetic energy and helicity, and will be solved in 
section V. 
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IVc. Quasi-normal approximation 

We now show that the rcsuhs in the above two subsections, can be also derived using a two-point-closure approxi- 
mation, also known as "quasi-normal" approximation. 

We start by observing that the autoinduction equation (|52|) . in Fourier space, reads: 

— ^7 — ^'^'^ W~]3 / (2 13 ^'3fc^jg'-^^(q) [£fers5n(k'-q)B„(k-k') - ersm^A;(k - k')B„(k' -q)] -r]k'^B^{k), 

(66) 
where Sijk is the totally antisymmetric tensor and summation over repeated indexes is understood. In quasi-normal 
approximation, the four-point magnetic correlator is decomposed, in terms of two-point correlator, as |23| : 

{B,{k)B,{k')Bk{ci)Bi{q'))^{B,{k)B,{k')){Bk{q)Bi{ci')) + {B,{k)Bk{q)}{B,{k')Bi^^^^^ 

(67) 
where (...) denotes an ensemble average. Multiplying Eq. ((66|) respectively by B*{k) and A*{k), and then averaging 
out, we get 



dt 
OHb 



-2?7efffc'£i3-Hai3fc'-HB, (68) 

-2r]^ffk^nB+4:aB£B- (69) 



dt 

As it is easy to check, the solution of the above equations in the case of magnetic fields with initial fractional helicity 
is given by Eqs. (J64p - (I65I) . so that the one-point- and two-point-closure approximations give exactly, and perhaps 
surprisingly, the same identical results. 

V. RESULTS 

In this section, we find the evolution laws of the magnetic intensity and correlation length in a radiation-dominated 
universe for the case of free-steaming neutrinos and photons, by solving Eqs. (|64 p - (j65p for some specific initial condi- 
tions. When working in a Minkowski spacetime, we take i = as the initial time for the sake of simplicity. 

Va. Initial conditions 

Magnetic spectrum- For the sake of simplicity, we assume that the initial magnetic energy spectrum is represented 
by the simple function 

fB(/c,0) = AsfcPexp(-2fc2£|), (70) 

where As and Ib are constants. For k -^ ig , the magnetic energy spectrum possesses a power-law behavior with 
exponent p which, to avoid infinities in the magnetic energy, we assume to be greater than 1. For the same reason 
(finiteness of energy), an exponential cut-off has been introduced. The quantity £b can be related to the initial 
correlation length, ^b(O), by Ib = r[(l +p)/2] CB(0)/[\/2r(p/2)], where T{x) is the Euler gamma function [56| . 

It is worth noting that the power-law behavior for the initial magnetic spectrum is predicted by many models 
of generation of cosmic magnetic field in the early universe. In particular, infiation-produced magnetic fields [lO[ 
usually exhibit a power-law spectrum when they re-enter the horizon. Also, most of the generating mechanisms after 
inflation [25] reposes on microphysical processes acting only on small scales. Therefore, the structure of the resulting 
magnetic field appears as a set of uncorrelated (Gaussian-distributed) eddies, which implies a fc^-like spectrum. To 
see this, suppose that the magnetic field at the initial time i = is a stochastic variable with gaussian distribution, 
7'[B(x,0)] ex exp[— (Stt/Ab) /d^a;B^(x, 0)] , where Xb is a constant. The magnetic energy spectrum, 

SBik, 0) = 2T:k^{\B(k, 0)n = 2nejd^xjd^y e^^^x-y) ^(i3,(x, 0) Sfc(y, 0)), (71) 

k—1 

is then proportional to /c^. In fact, observing that 

/Rl n^T, / m\ /P[^fe(^. 0)]^[B(z, 0) ] i3.(x, 0)B,{y, 0) A^ ^ _,, . ,_ 

(S.(x,0)B,(y,0)) = Jv[B,{zMnB{z,0)] = ^ ^.. '^(x - y), (72) 
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where in the last equality we used the properties of Gaussian integral [57|, we get Ssik, 0) = XBk"^ ■ 

Recently enough [58|, however, it has been shown that the analyticity of the magnetic field correlator Bi(jc)Bj{y) 



defined on a compact support, together with the divergenceless of B, implies that the spectral index p has to be even 
and equal or larger than 4. Moreover, in [5C|, it has been shown that, starting from suitable initial conditions which 
should reflect the production mechanism of a magnetic field from bubble collision in a first-order phase transition, 
a Batchelor spectrum £b oc fc'' is established after a short time interval and at small wavenumbers. Nevertheless, 
we can always assume that on very large scales (very small wavenumbers k) the initial spectrum satisfies the above 
requirements, but that on smaller scales, where the gross of the magnetic energy is stored, it behaves like k^. In 
any case, in the following we leave the index p as a free parameter, privileging the cases p — 2,3,4 when we show 
graphically our results (for the choice p = 3, see below). 

Drag coefficient.- We assume that the drag coefficient a scales in time following the simple law 

a(i)=a(0)[l+r/7(0r, (73) 

where r is the normalized time 

Teddy yyj) 

and 7(0) and a are constant. The above parametrization is useful since a{t) ~ a(Q) for t <C T^ddyi^), when the system 
is, as we will see below, in a quiescent phase [i.e. the integral quantities like EB{t) and S,B{t) remain almost constant in 
time], and a{t) follows asymptotically a power-law behavior Q;(0)[r/7(0)]° for large times, when magnetohydrodynamic 
effects operate in changing the state of the system. Using the results in [53i . |54] , we find that in comoving variables 
a = —3 for the case of free-steaming photons, while a = —4 for the case of free-steaming neutrinos. The constant 7(0) 
is generally different from unity, and its explicit expression will be derived in section Vh. 

Setup- Numerical integration of dragged MHD equations has been performed only in |53l l54l| . There, the simple 
case of constant dragged coefficient was analyzed. Also, one of the cases discussed in 53|, [SJ], was that with p = i. 
Hence, in the light of the above discussions, we consider the three cases 

{(2,-3), Gaussian, photon case, 
(3,0), Banerjee-Jedamzik case, (75) 

(4, —4), causal, neutrino case, 

for the non-helical case, and the three cases 

{(2, —3, 10^^), Gaussian, photon case, 
(3,0,1), Banerjee-Jedamzik case, (76) 

(4, —4, 10^^^), causal, neutrino case, 

for the helical case. Note that in the helical, Banerjee-Jedamzik case, the magnetic field is taken to be maximally 
helical, as it is in the numerical analysis of [53|, l5J|. Also, the very small value h = 10""'^^ is taken to clearly show the 
transition from the selective decay phase to the inverse cascade phase to be discussed in section Vd. 

Vb. Master equations 

Inserting Eqs. (|64| and (|65|) in Eqs. (O and ^ we find, respectively, 

EB{t)=EBmi + CLs)-'^ l<^l+phB^^=k==^2], (77) 

HBit) = HBmi + CLs)-^ I *3 + 7^ ;?- ^=^ *4 ) , (78) 

with 
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where iFi{a,b;z) is the Kummer confluent hypergeometric function [56|, and we have introduced the quantities 
Ca(i) = C(t)/^B, Cdiss(t) = edissit)/£B, and Cb = ^5(0)7^5. Also, inserting Eq. dH in Eq. ^ we get 



?B(t) = eB(o)(i + ciss) 



.2 Nl/2 '^^ \/l+C 



It is useful, for the following discussion, to define accurately the magnetic Reynolds number and the eddy turnover 
time as 

-Q WrmsSiS t,B /oo\ 

-n-e_B = , 7'cddy — , (O^) 

where we used the magnetic correlation length and the root-mean-square value of the velocity field, v'^^^ — 
Jd^xv^{x,t) = 2Ey, as typical length scale and velocity, respectively. With the aid of the above definitions, the 
integro-differential equations ([77]) and ([751) for Eb and Hb can be transformed into ordinary differential equations for 
the quantities Cdiss and Cq. We have 

^._k„o„„R,,(„)[i + ./,(„,!-• ^|f|, (84) 

with Eb and Hb as a function of ^diss and Ca given by Eqs. (1771) and (j78p . respectively. The constant 

m = ^1^ (85) 

Tcddy(0)r(0) 

is of order unity [see Eq. ([T9| ]. In section Vf, we will found that this order-unity constant is indeed a parameter which 
weakly depends on p. 

Equations ((77)) - ((75)) and ((55)) - ((M| are our master equations in the study of the evolution of the magnetic field. 

Vc. Solutions: non-helical case 

In the non-helical case, Hb — 0, the solution of Eqs. (|83| - (|84p is such that Ca{t) = 0, that is HB{t) — for all times. 
Then, from Eqs. dTT]), dH]), and dUl) we get 

EB{t)=EB{Q){^ + CLsr'^, (86) 

ut)^uo)a + cLs)'^', (87) 

£B{k,t) = £B{k,0)exp{-2eelClJ. (88) 

For large magnetic Reynolds numbers, the first term in the right-hand-side of Eq. (|55)) can be neglected with respect 
to the second one. In this case, the solution of Eq. (|55| is 

CLsir) = {co{[l + r/jiO)]'-^ - 1} + 1}^ - 1, (89) 

where cq is given in Appendix A. For r ^ 1, we have 

Cdiss(^)-ClT^fe^, (90) 



where ci is given in Appendix A, so that 



Eb{t) ~ C2Eb{0)t~^''^^+p^"\ (91) 

CB(r)~ 03^5(0) rfet, (92) 

£B{k,t) ~£:B(fc,0)exp[-C4fc^i^^j, (93) 
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where 02,03,04 are given in Appendix A. Equations ([9T|) . ([92t . and (l93l) . are in agreement with Eqs. (|45)) . (|47l) . and 
(|43)) . respectively. We deduce that the mean- field-theory approach and the simple scaling arguments applied to MHD 
equations are compatible. The mean-field-theory approach gives us two more pieces of information with respect to 
scaling arguments. Firstly, the time when the system enters into the scaling regime is, looking at Eqs. (1551) . (|87p . 
(|55)) . and (|5^. approximatively given by the initial eddy turnover time times 7(0). Secondly, comparing Eq. (1^1) with 
Eq. ((321) we get the expression for the scaling function ij^Bix): 

^B{x)^e-''^\ (94) 

In Fig. 1, we plot the spectrum of the magnetic energy at different times for the three cases in Eq. (175]) [in all 
figures of this paper, we take 7(0) = 1]. It is clear that the decay of the magnetic energy and the growth of the 
magnetic correlation length (which is, roughly speaking, the inverse of the wavenumber corresponding to the peak in 



the magnetic spectrum) advance through selective decay, a phenomenon well described in 38| in the case of turbulent 
MHD (see also [23|). There is no direct transfer of magnetic energy from small scales to large scales but modes with 
larger wavenumbers just decay faster than those whose wavenumbers are small. Accordingly, as the time goes on, the 
magnetic field survives only on larger and larger scales, as it is evident in Fig. 1. 

The dotted lines in Fig 1 represent the magnetic energy spectra at a certain time i/Tcddy(0) = Tg (to be defined in 
section Vd) when, if the field were helical, the magnetic field would enter in a phase of inverse cascade (see below) . 
In the case at hand, the magnetic field has no helicity and, indeed, nothing special happens when r reaches ts- 

Vd. Solutions: helical case 

In the helical case, the dynamo and dissipation lengths, as well as the ratio Ca/Cdiss, are increasing function of time. 
This is inferred by the results of numerical integration of Eqs. ([83)) - ([84)) as shown, for example, in the left panel of 
Fig. 2. [Here, and in the following, we consider the case of a very high magnetic Reynolds number, Res = 10^^]. In 
the Appendix B, we show that there are two different regimes depending on the value of Ca/Cdiss- For Ca/Cdiss «C 1 the 
system behaves as if the magnetic helicity were zero, so that the system evolves by selective decay. The asymptotic 
solutions are then the same as obtained in section Vc. For Ca/Cdiss ^ 1, instead, the mechanism operating in the 
evolution of the system is such that energy is transferred from small to large scales. In turbulent MHD, this mechanism 
is known as inverse cascade and will be discussed, in more detail, later on in this subsection. 

The following asymptotic (t — > 00) solutions of Eqs. ([83]) and (|84|) are derived in the Appendix B in the case of 
very large magnetic Reynolds numbers: 

Cdiss(r)~C5(lnr)i/V(i-'^)/3, (95) 

Ca(r)^C6(lnr)2/3r(i-)/3, (96) 

and 

Eb (t) ~ C7 Eb (0) (In r)i/3 t-'^^-)/^ , (97) 

CB(r)^C8CB(0)(lnT)-i/V(i-")/3, (98) 

where 05,06,07,03 are given in Appendix A. The above relations have been obtained supposing that the magnetic 
helicity is almost constant during the evolution of the system. This assumption is well justified by the results of the 
numerical integration of Eqs. (|83l) and (|M1) . as shown, for example, in the right panel of Fig. 2. 

In the upper panels of Fig. 3, we show the magnetic energy and correlation length as a function of time for 
the three cases in Eq. (l76l) . Dotted lines are the asymptotic expansions (|97| and (|98)) for t > te and t > t^, 
respectively, and the analytical results (|86l) and (|87l) [with Cdiss given by Eq. (|89|) ] for < r < r^ and < r < t^, 
respectively. Here, te and r^ are the times when, approximatively, the integral quantities Eb and ^b enter into the 
inverse cascade regime. We can find them by simply matching the asymptotic solutions in the two different regimes, 
namely equating Eqs. dH]) and ^ to Eqs. ^ and dMl), respectively. We find, te ~ (07/02)"' [ln(o7/o2)"«]~'^^ 

and Tj ~ (os/os)"^' [ln(o8/o3)^^'] '' , where q ~ 3(3 +p)/[2p{l — a)]. In the case Hb ^ 1, the above results can be 
approximated as 

rEO,c,{\nh~^y^'h-^"^/', (99) 

T^^c^o{\nh^')''^'h~^"^/', (100) 
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FIG. 1: Energy spectra in the non-helical case for (p, a) — (2,-3) (upper left panel), (p, a) = (3,0) (upper right panel) and 
{p,a) — (4, —4) (lower panel). Dashed lines correspond to the initial spectra. Continuous lines correspond, from left to right, 
to t/reddy(0) = 1, 10, 30, 10^ 3 X 10^ lO'*, 3 x lO" (upper left panel), f/roddy(0) = 1, 10, 10^ 10^, lO", 10^ 10®, lO'^ (upper right 
panel), and i/rcddy(0) = 1,10,10^,10^,3 x 10'*,10'',3 x 10*^ (lower panel). Dotted lines correspond to i/rcddy(0) ~ 2.3 x 10^ 
(upper left panel), and i/roddy(0) ~ 3.5 x 10^ (lower panel). 



where cg,cio are given in Appendix A. 

We also note that, in the inverse cascade regime, the product of magnetic energy and correlation length is a 
quasi-conserved quantity. In fact, as we show in the Appendix B, we have 



EBir)^Bir) 



Hsir) 



(101) 



From the above equation, we deduce that a magnetic field with fractional helicity becomes maximally helical approx- 
imatively after the system enters into the inverse cascade phase. The time when this happens, th, can be found by 
matching the product of asymptotic solutions Eqs. (pT|) -(|92 |) and (|101[) . We get 



TH ^ Cllh 



'2g/3 



(102) 



where cn is given in Appendix A. In the lower panel of Fig. 3, we show the quantity 2EbS,b/ Hb as a function of time 
for the tree cases in Eq. (|76)) . The asymptotic expansions for t < th and t > th are practically indistinguishable 
from the numerical solutions. 

In the helical case, the evolution laws of magnetic energy and correlation length do not depend on the index of 
the initial magnetic energy spectrum. Moreover, apart logarithmic corrections, and for the case of constant drag 
coefficient (o = 0), we find the t~^'^ and t^'^ laws which are indeed observed in numerical simulations of dragged 
MHD [5i,[54|. 
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FIG. 2: Evolution of the a-length normalized to the dissipation length (left panel), and magnetic helicity (right panel), for the 
tree cases in Eq. (|76l) . 
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(indistinguishable from the continuous lines in the lower panel). 



It is evident from all the above figures that the (approximate) analytical expansions fit very well the numerical 
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solutions. Moreover, the magnetic helicity is, as expected, an almost constant function of time. Because of the 
conservation of the magnetic helicity, small-scale modes are not dissipated during the decay but their energy is 
transferred to larger scales. This process of inverse cascade is manifest in Fig. 4 (see, in particular, the magnetic 
spectra). It is worth noting that, as explained above, this mechanism begins to operates only after a certain time, 
that is when the selective decay mechanism ends. We can suitably define the time when the inverse cascade begins, 
T£ , as the time when the maximum of the magnetic spectrum meets the initial spectrum (see the magnetic spectra in 
Fig. 4). Assuming for simplicity /i^ ^ 1, we get (see Appendix C) 

T£ ~ Ci2/1b''/' (103) 

where ci2 is given in Appendix A. Looking at Eqs. (|74)) . ([99]), (jlOOp . (|102p . (|103p . and roughly speaking, we can say 
that the time when system enters into the inverse cascade regime is iinvcas ^ hg Piiypy ° Teddy (0). Finally, we 
note that magnetohydrodynamic effects do not change the characteristics of the magnetic field on scale well below 
the integral scale (the magnetic correlation length) since small-scale modes need longer times to be processed by the 
inverse cascade. This is true for all the cases we have analyzed, although it is visible just in the upper left panel of 
Fig. 4, where it is evident that the magnetic energy spectrum retains its initial power law form at very large scales 
(in the other cases, the magnetic energy and helicity spectra are not shown for very small wavenumbers) . 

Ve. Solutions: maximally helical case 

In the case where the initial magnetic field is maximally helical, ft, = 1, we get two exact analytical results. 
Firstly, from Eqs. ([M)) and (p5|) we find that the magnetic field remains maximally helical for all times, HB{k,t) = 
2k~^£B{k,t). Secondly, integrating the above expression in k we have EB{t)^B{t) — HB{t)/2, a result usually used 
in the literature for the evolution of a (maximally) helical magnetic field, but never fully justified. 

Vf. Kinetic energy 

Starting from Eq. (1201) and taking into account Eq. (|73p . it is straightforward to find the asymptotic expansions of 
Ey for T > 1: 

„ . > \ T 3+p non-helical, /^r^A\ 

^^^^^ "^ \ n M/3 -^<H±^ ;. r 1 (^"^) 

[ (In t)^/'^t 3 , helical, 

where we used Eq. (|91]) for the non-helical case, and Eq. (|97l) for the helical case. Equation ()104p is in agreement 
withEq. (|46| . 

Moreover, evaluating Eq. ([20)1 at the time r = 0, and taking into account Eqs. (IT7|) . ([55]) . (1551) . and Eq. (I55|) in the 
limit of large magnetic Reynolds number, we straightforwardly find that Ey{Q) — [(1 +p)/12] £^B{0)^S{0)'^Ey{0), from 
which we get 

12 V/^ 1 



-5(0)- -— — . (105) 



Finally, starting from the definition of (5(0) in Eq. ([55]) . and using Eqs. ((T7)) and ([5^ evaluated at r = 0, we find 

^b(O) 



EM 



L<5(0)«(0)^b(0) 



(106) 



which relates the kinetic energy at the onset of the dragged phase to the initial values of the drag coefficient, the 
magnetic field energy, and the correlation length. 

Vg. Exiting the dragged phase 

Neutrino and photon mean free paths- The dragged phase is characterized by the condition that the comoving 
mean free path of neutrinos, t^f„, and/or photons, i^L, is much greater than the comoving magnetic correlation 
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FIG. 4: Upper panels. Magnetic energy (left panel) and helicity (right panel) spectra in the helical case for (p, a, h) = 
(2, —3, 10~^). Dashed lines correspond to the initial spectra, while continuous lines correspond, from left to right, to f/reddy(0) = 
1, 10,30, 10^3 X 10^, 10^,3 X 10*. Dotted lines correspond to i/reddy(0) = r£ ~ 2.3 x 10^. Middle panels. As in the upper 
panels but for (p, a, h) — (3, 0, 1) and t/r(,ddy(0) = 1, 10, 10^, 10^, 10*, 10^, 10^, 10^. Lower panels. As in the upper panels but 
for (p, o, h) = (4, -4, 10^-^^), f/roddy(0) = 1, 10, 10^ 10^ 3 x 10*, 10^ 3 x 10^ and re ~ 3.5 x 10^ 



length. Therefore, it is not obvious that this condition in maintained during the evolution of the system, since ^b 
and both £j^/p and t2l^ evolve in time. Observing that t^^^ oc T^'^ (x R^ (x ¥ and ^J^jp ex T^^ oc R"^ (x P in 
radiation-dominated era 53l |54|, we have £^L ex t"* and ^^1 oc r^. Here, we used the results in section lid, the 
fact that R ex <^" in radiation era (with t being the cosmic time), and the fact that the temperature T and the 
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scale factor R are related by i? ex g~g' T~^ [595, '^ith g^,s being the effective number of entropy degrees of freedom 
at the temperature T. On the other hand, the maximum growth for the correlation length happens for the helical 
case and for the neutrino free-streaming case (a = —4) [see Eqs. (|92|) and (|98l) ]. Namely, ^b increases at most as 
^B oc (Inr)"^/'^ T^/'^. As a consequence, if the system is in a dragged phase such that ^b ^ ^mfp' ^mfp' ^^ ^^^^ remain 
in it for all times. 

However, the evolution laws for the magnetic energy and correlation length has been obtained in the equilibrium 
state, where the kinetic Reynolds number is much smaller than unity. Re <^ 1. This quantity generally evolves 
in time as the typical velocity and correlation length of the fluid motion, and the drag parameter change in time. 
Consequently, there could exist a time when the condition Re ^ 1 ceases to be satisfied and the system exits the 
dragged phase. 

Kinetic Reynolds number.- Defining the kinetic Reynolds number accurately as 

Re = ^ , (107) 

we find Re = {—EbY^^ ^^^ a^^^^ , where we used Eq. (pil)) . The asymptotic expansions (for r ^ 1) of the kinetic 
Reynolds number in the non-helical and helical cases, and for large magnetic Reynolds numbers, are 

■D / \ / T"(i+''\ non-helical, ,_,_^. 

^^(")°^|(lnr)V2r-a+'^), helical, (^^^^ 

where we used the fact that ^v{t) oc ^B{t) [see the first relation in Eq. (|47l) ]. Eq. ([73l) . Eqs. (|91]) -([92 |) for the non-helical 
case, and Eqs. (|97| -([98 |) for the helical case. From Eq. ()108p we see that in the physical cases of interest, a = —3, —4, 
the kinetic Reynolds number is an increasing function of time. Therefore, there will be a time tcxit when the system 
leaves the dragged phase characterized by Re <C 1. This time can be operatively defined by the condition Re(texit) — 1- 
Finally, let us remember that the kinetic Reynolds number is approximately equal to the F ratio, Re(t) ^ F(i) [see 
Eq. (HH)], so that Eq. (fTOS]) is in agreement with Eq. (pS]) . 



Vh. Evolution laws in expanding universe 

Returning to the "tilde" notation and indicating the time when the dragged phase begins as ti , equation ([22l) reads 

a(i) = aiU) fl) . (109) 

Introducing the new variable f as 

r = ^T-, (110) 

Toddy [ti ) 

we can write Eq. (|109l) as a{t) = a{ti)[T/j(ti)]"' with 7(ii) = ii / Teddy (^0- ^^ radiation-dominated era, it is easy to see 
that j{ti) can be expressed as 7(ii) = [c!_fr(ii)/?B, phys(^i)] i'rms(^i): where dn = 2i is the length of the Hubble horizon 
at the time t, with t being now the cosmic time. Therefore, in the notation of section Va, we have 7(0) = NiVi, where 
we have introduced the initial number of magnetic domains in the horizon 

t,B,phys[(j) 

and we have defined Vi — Wrms(O) for the sake of simplicity. Moreover, the normalized time f in the radiation-dominated 
era is easily found to be 

f = N,v,^. (112) 



Finally, using Eq. (|106p . we can relate Vi to the initial values of the magnetic field strength, B^ = B 
Jd^x'B'^{x,t) — 2Eb, the physical magnetic correlation length, ^s^phys — R^b, and drag coefficient, as 



rins 



Oi(-rC.i) 1;b, phys(^i) 





R. 


<R< 


Ri 




Ri 


<R< 


i?2 


(l-a)/6 


R 


>R2, 
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where kq is given in Appendix A. [We note that in the notation in which the quantity p + P is not taken to be 1, the 
magnetic field strength B{Ri) in Eq. (|113p should be replaced by B{Ri)/ ^/p{Ri)'^~'P{Ri)] 

Non-helical case- Translating the results previously found in the case of a Minkowski spacetime to the case of a 
flat Friedmann universe, we get 

R'BjR) ,,[^' <_,,,, R^<R<Ru 

for the comoving magnetic field strength, and 

URT) ^ j «2 {N,v.y/<-'^p^ (f )^ ,R>Ri, ^^^^^ 

for the comoving magnetic correlation length. Here, 

^ c, KsiN^v^y'/^'-'^K (116) 

Hi 

and Ki, K2, H3 are given in Appendix A. In Eq. (|116p . we have assumed that K3{NiVi)~^^^^~'^'' > 1. If this is not the 
case, Ri is (approximatively) equal to Ri. 

Helical case- In the case of helical magnetic fields, we find 

Hg^ J MM.r-^ (f )-^ «.<«<«. ,UT, 

for the comoving magnetic field strength, and 

{1, Ri ^ R ^ Ri, 

«2 {N,v,y/(^+p) U)^ , Ri<R<R'2, (118) 

1 /, / \-l/3 / \(l-a)/3 

«;5/iB^'(A^.«.)^/'(lnA^,«,f) (f) ,R>R'2, 

for the comoving magnetic correlation length. Here, 

R2 ^6 /, \-q/3 R2 ^7 n \2n/3 /'11n^ 

— ~-^w(lnw «/^ — ~ -^cj Inw) «/^ (119) 

Ri NiVi Ri NiVi 

with uj — [NiVi) " " /i^ "^ , and K4, Ats, Kg, ^7 are given in Appendix A. In the above equations, we have assumed 
that uj is much greater than unity. 

VI. DRAGGED MHD IN MATTER-DOMINATED UNIVERSE: FREE-STREAMING PHOTONS AND 

HYDROGEN ATOMS 

We now discuss the evolution of a magnetic field in matter-dominated universe, when the drag term in the MHD 
equations is ruled by the free-streaming properties of photons or hydrogen atoms. 

Via. Evolution laws in supercomoving variables 

The drag coefficients for photons and hydrogen atoms evolve as a function of the temperature as 53|, [54 1 a-y (T) ex T^ 
and aniT) ex T"^, respectively. In supercomoving variables (namely in Minkowski spacetime), their evolution with 
time is then straightforwardly given by 

a(f)=a(0)e'^^/T(°), (120) 
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where, taking into account the results of section lid, we used the fact that the normahzed time f in Eq. (IllOp is 



1 .. , ^ 

— JNiVi In — 

2 R, 



in matter-dominated era, when the scale factor evolves as i? ex t^^^ . In Eq. (|120p . 

—3, photons, 
—5, hydrogens, 



(121) 



(122) 



and 7(0) = NiVi/3. The initial number of magnetic domains per horizon, Ni, is the same as in Eq. (Illip . the length 
of the Hubble horizon in matter-dominated era being dn — 3t (with t the cosmic time), and Vi — fims(O) is given, as 
in the case of radiation-dominated universe, by Eq. (jllSp . 

In supercomoving variables, the master equation (j83p reads, for the case at hand, 



"Cdiss 






^ (5(0) e""^/^(") ^^^^^ 



df Rei3(0) 3 ' ' EsiO)' 

Defining the new variable 7(0)[e~°'^/''(°) ^ 1] = NiVi{R/ RiY/"^ /i, the above equation becomes 



dClss 



a 



df Res(O) 3 



^(5(0)[l+f/7(0)]-('^+i)^^^ 



Eb{0) 



(123) 



(124) 



Comparing Eq. (I124p with Eq. (|83)) , we see that the evolution of the magnetic energy and correlation length can be 
derived by applying the substitutions r — > f and a — > a -|- 1 in the solutions found in the case of a Minkowski spacetime 
in sections Vc, Vd, and Ve. 



VIb. Evolution laws in expanding universe 

Non-helical case- Using the results in the above subsection, we straightforwardly find the evolution laws for the 
the comoving magnetic field strength. 



R^B{R) 



1, 



R^B{R,)- yn,{N,v,)-^^ (fy 
and for the comoving magnetic correlation length, 

^b{R) „ 



3o(l + p) 
3 + p) 



^B{R^) 



Kg 



where 



(iv^„^)i/(3+p) (^y 



R\ ,-^ , _2_ 

— - ~ Kl0(iViWi)3a^ 



R^<R<Rl, 

R>Ri, 



Ri ^ R ^ Rly 

R>Ri, 



(125) 



(126) 



(127) 



and Kg, '^9: '^lo are given in Appendix A. In Eq. (|127p . we have assumed that Kio(A'iWi)^/'^'^°) > 1. If this is not the 
case, then i?i ~ Ri. 

Helical case.- In the helical case, we get 



R^BjR) 
RlB{R) 



1, Ri ^ R ^ Ri, 

3a(l+p) 

nsiN,v,r-^{^y'' , Ri<R<R2, 

, ,„ r / s 3/2"|l/^ / sa/4 



(128) 
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for the comoving magnetic field strength, and 

' 1, 



«9 (iV,U,)l/(^+P) (f ) 



'2(3+p) 



Hi ^ R ^ Ri , 
Ri ^R^R2, 



^BiRi) ~ 
for the comoving magnetic correlation length, where 



IniTV.t;, f 



3/2' 



-1/3 



-a/2 



(129) 



[-i-j ^ R^R'2: 



Ra 
Ri 



«13 



a>(lncj)-^/3 



2/3 






«14 



(I>(lncj)2'?/3 



2/3 



(130) 



with oj — [NiViY^^^i ' ""hg '^ /3 and q — — 3(3+p)/(2pa), and kh, K12, K13, K14 are given in Appendix A. In the above 
equations, we have assumed that w 3> 1. 

Vic. Exiting the dragged phase 



mfp 



Photon mean free path- In the matter-dominated era, the comoving mean free path of photons scales in time as 
ex T^^ oc R^ . On the other hand, the maximum growth for the correlation length happens for helical fields in 
the case of free-streaming photons. Accordingly, ^b increases at most as £,b oc (lni?'^/^)^^/'^i?^. As a consequence, if 
the system is in a dragged phase such that ^b ^ ^mfp' ^* ^^^^ remain in it for all times. 

Hydrogen atoms mean free path- The comoving mean free path of hydrogen atoms evolves as r^l oc T~^ ex R^. The 
correlation length in the hydrogen atoms free-streaming case, instead, increases as £,b oc i?i5/P(3+?')l in the non-helical 
case, and £,b oc (lni?'^/^)~^/^i?^/^ in the helical case. As a consequence, the system will eventually exit the dragged 
phase in the helical case, and in the non-helical case for the case p < 4/3 (we remember that we are assuming p > 1). 

Let us call ti^exit the time when this happens, which can be defined by the condition ^^(ti^exit) — ^^fp(ii,exit)- 

Kinetic Reynolds number- The asymptotic expansions of the kinetic Reynolds number in the non- helical and helical 
cases, and for large magnetic Reynolds numbers, are 



r ^-3a/2 
MR)'^{ A__^3/2^1/2 3a/2 



Ini?^ 



non-helical, 
helical. 



(131) 



From the above equation, we sec that the kinetic Reynolds number is an increasing function of time in the physical 
cases of interest, a = —3, —5. Therefore, there will be a time i2,cxit when the system will leave the dragged phase 
characterized by Re <^ 1. This time can be operatively defined by the condition Re(i2,exit) = 1- Accordingly, and in 
the cases of physical interest p — 2,A, we have that the dragged phase terminates at the time toxit given by 



i2,Gxit J photon case, 

i2,oxit , hydrogen, non-helical case, 

min{ti_oxit , t2,cxit}, hydrogen, hehcal case. 

Finally, we observe that, as in the radiation-dominated era, we find that Re ^ F, in agreement with Eq. (j25p . 



(132) 



VII. JUSTIFYING THE BANERJEE AND JEDAMZIK'S SCALING ARGUMENTS 



In this section, we reconsider some scaling arguments used by Banerjee and Jedamzik in 53|, [5^ to explain their 
results obtained by direct numerical integration of dragged MHD equations. Our goal here, it to show that some 
working hypotheses used in \53i, iM] find indeed a justification in the light of our analytical results obtained in sections 
II, V and VI. 

Master^uation- In order to discuss the Banerjee and Jedamzik's scaling arguments, let us introduce the notations 
used in |53|, |54|. We define the magnetic energy in the mode k, EB{k,t), by EB{t) = j^dkk~^EB{k,t). This is 



related to the magnetic energy spectrum by EB{k,t) = kSsik^t). Accordingly, we can define the magnetic energy 
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on the scale I ~ 1/fc as EB{l,t) — {l/l)£B{^/l,t)- At the initial time t — 0, the magnetic energy on the scale I is 



proportional to a power of /, Eb{1, 0) oc Z ", where n = 1 + p in the notation of [53|, |54| 



In [53J,l5J], it is assumed (without a full justification) that the evolution of the magnetic energy on the integral scale 

(133) 



I sit) (namely the scale under which the magnetic energy is dissipated) is ruled by the simple scaling equation 

dEsilB) Eb{Ib) 



dt 



'eddy 



where Eb{Ib) — EB{lB{t),t). The above equation can be justified in the light of our results in section lib. In fact, 
from Eq. (gOl), it follows that 

^^-2^^, (134) 

dt Teddy 

where EB{t) is the magnetic energy. In order to obtain the above equation, we eliminated Ey in Eq. (j20p in favor 
of r and then used the fact that aT ~ l/roddy [see Eq. (TTOl) ]. We now use the fact that the magnetic spectrum, 
in both the helical and non-helical cases, is peaked at a wavenumber fcniax(0 ^ l/^s(i) ^ ^/Isit), as it is not 
hard to see using the expressions for the magnetic energy spectrum and correlation length in section V, so that 
EB{t) ~ J^dk£B{k,t) ^ fcmax£B(fcmax,0 ~ Eb{Ib)- This last observation, when combined with Eq. (|134p . gives 

Eq. ((I33D. 

Now, we can transform Eq. (jl33p into an equation which does not involve the bulk field v (which enters in the 
equation through the quantity Teddy)- To this end, we use Eq. (j38p and approximate the nabla operator by V ~ 1/?b- 
We find that the typical velocity v and the typical magnetic field strength B are related through v ^ B"^ /{uIb) ^ 
Eb/{olIb)- In turn, this gives for the eddy turnover time Teddy ^ ^%/{o^B). Inserting the above equation in Eq. (jl33p 
we get 

dEBJls) EIjIb) ,.„,, 

This is the master equation used in [53l l54l | to derive the evolution laws for the magnetic energy and correlation 
length. As we have just shown, this equation is a direct consequence of Eqs. (dU and (1^ . 

Non-helical case- In [53J, |54|, it is assumed that the magnetic energy on the integral scale Ib is related to the 



integral scale itself by a simple power law 

Eb{Ib) oc ?b". (136) 

This is indeed legitimate by the fact that, as we have shown in sections V and VI, the evolution of the magnetic 
spectrum in the non-helical case proceeds through selective decay, so that the initial magnetic spectrum retains its 
form (for length greater than the integral scale) for all times. 

From Eq. (fT36l) . we have Ib oc EBihy^^" that, when inserted in Eq. (fT35l) . gives ^^(t) ex [/dt/a(i)] ""^^^^"^ 
and, accordingly, lB{t) oc [J dt/a{t)\ . Remembering that the index n is related to our index p through 

n = 1 + p, and taking into account the expressions for a as a function of time, we get EB{t) ex i~(i~'*)(i+p)/(3+p) 
and £,B(t) ex f(i~°)/(3+p) jjj ii^q case of free-streaming neutrinos and photons in radiation-dominated universe, and 
EB{t) ex e'^'^^+P-'/'^'^+P^l * and iB{t) ex e~[''/^'^+P^l * in the case of free-streaming photons and hydrogen atoms in matter- 
dominated universe. 

The above evolution laws are in agreement with those derived in mean- field approximation (see sections V and VI) . 

Helical case- In the helical case, the magnetic helicity is related to the magnetic energy on the integral scale Ib 
through Eb(Ib)Ib — Hb{Ib)- Since the magnetic helicity on the integral scale Ib is a quasi-conserved quantity, 
Hb{Ib) — constant, we have Ib oc Eb{Ib)~^- Inserting the previous equation in Eq. (|135p . we obtain EB^t) ex 
[/(ii/a(t)] and, consequently, lB{t) oc \^J dt/a{t)\ . Writing a as a function of time, we get EB{t) ex t^(i~'^)/3 

and ^B(i) oc i(i~'^)/3 for free-streaming neutrinos and photons in radiation-dominated era, and EB{t) ex e*^"/'^-'* and 
S,B{t) oc e~^""'* for free-streaming photons and hydrogen atoms in matter-dominated era. 

Also in the helical case, the above evolution laws are in agreement with those derived in mean- field approximation 
(see sections V and VI). 
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VIII. FREELY-DECAYING TURBULENT MAGNETIC FIELDS 

The problem of determining the evolution properties of freely-decaying magnetic fields in MHD turbulence, in 
contrast to the case of dragged MHD, has been deeply and widely discussed in the literature. Essentially, there are 
been five ways to tackle it: i) Applying scaling arguments to turbulent MHD equations 3l|, l35|, IST], |43( ; ii) Using 



simple toy models to "mimic" the full turbulent MHD equations [34, 36;, 43]; in) Solving the turbulent MHD equations 
in Eddy-Damped Quasi-Normal Markovian (EDQNM) approximation [38|; iv) Studying MHD turbulence equations 



in mean- field-theory approximation 48 1; v) Solving numerically the turbulent MHD equations [39l. kol. l43l 1471 149|- 



The invariance of turbulent MHD equations under scaling transformations has been firstly used by Olesen [35| to 
study the evolution properties of the magnetic spectrum and correlation length. Although this is a very simple and 
elegant way to study MHD turbulence, this approach lives some unresolved questions, such as when the system enters 
into the scaling regime, what is exactly the shape of magnetic spectrum, etc. 

The most famous toy model for hydrodynamic turbulence is the so-called shell model based on the original idea by 
Desnyanski and Novikov [60|, developed later by Glcdzcr 6JJ, Ohkitani and Yamada [621, ^''^'^ ^^'^ know as GOY- 
model. Subsequently, it was generalized by Gloaguen et al. to include the case of magnetic turbulence [63|. The 
basic idea on which repose the GOY-model is that the interaction due to non-linear terms in the MHD equations are 
local in /c-space. Then, convolutions in MHD equations are approximated by sums over the nearest and next nearest 
neighbors in a discretized wavenumber space (the shell space). Moreover, the shell-model equations mimic the full 
MHD equations retaining their main characteristic, such as energy conservation in the case of null dissipation, phase 
space conservation, etc. Also, there exists a continuous version of discretized shell model in which the distance between 
shells goes to zero. It was firstly introduced by Parisi [GJ] in a hydrodynamical context and later extended to the case of 
MHD turbulence by Brandenburg, Enqvist, and Olesen [36J. Although many properties of turbulence, such as energy 
transfer, general spectral properties, etc., have been studied successfully using shell models, the physical validity of 
such toy models remains an open question (for a review on shell models of magnetohydrodynamic turbulence, see [65|). 

The EDQNM approximation is the most widely two-point-closure approximation used for studying hydrodynamical 
turbulence. It was first proposed by Millionshtchikov [66], developed by Proudman and Reid [67|, and extended to 
MHD turbulence by Pouquet et al. [68|. It consists in replacing correlation functions by suitable statistical aver- 
ages [44] . neglecting (conveniently) high-order correlators, and using a phenomenological expressions for the so-called 
"eddy damping rate" . Even if the resulting equations are more simple than the full MHD equations, they can be 
solved only numerically and, moreover, some unresolved questions concerning EDQNM approximation still exist [23|. 

In [48^, analytical results for the evolution of the magnetic energy density and correlation length were obtained 
in mean-field-theory approximation. Below, we briefiy review such an approach to the study of turbulent MHD 
equations, and extend the results obtained in non-expanding universe to the case of a flat Friedmann spacetime. 

A direct integration of the full set of turbulent MHD equations is certainly the best way to study the dynamics of 
freely-decaying MHD turbulence. However, due to the complexity of such highly non-linear equations, the numerical 
results in the literature are often contrasting 

In Tab. 1, indeed, we show some results found in the literature concerning the evolution laws of magnetic energy 
and correlation length in freely-decaying MHD turbulence in Minkowski spacetime. In all cases, the assumed initial 
magnetic energy spectrum has a simple power-law behavior, and the dissipation parameters are constants. Obviously, 
one should keep in mind that these results were obtained by using different initial conditions, such as initial magnetic 
Reynolds number, initial F-ratio, an so on. As it is evident from the Table, the situation is still unclear. Nevertheless, 
we see that convergence in the results, at least for the most recent works and for the case of helical magnetic fields 
seems achieved (see, in particular, Refs. |48ll5Cll5lU53ll54|). 



Villa. Evolution laws in non-expanding universe 

Prior to recombination, namely for temperatures T > Tj-cc — 0.3 eV [59|, the equations of turbulent MHD are given 
by Eqs. ((T]) and ^ by the replacement 

- av -> i^V^v, (137) 

where ly is the kinematic viscosity. After recombination, instead, and in the tight-coupling regime between ions and 
neutrals particle species, the turbulent MHD equations are given by Eqs. H]) and ^ by the replacement 53ll54| 



av^ Vx(FlxB), (138) 

Piain 
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TABLE I: Evolution laws for the magnetic energy, Esit) oc f^*^, and magnetic correlation length, i^s(i) oc i*", in helical and 
non-helical freely-decaying niagnetohydrodynamic turbulence. All results were obtained assuming an initial power-law for the 
magnetic energy spectrum, £'_b(A:,0) oc fc*", and constant dissipation parameters. The asterisk refers to analytical solutions of 
continuous shell model. 



Reference 


Method 


P 


O" 


<; 


Helical 


[35] 


Scaling 


P 


- 


2/(3 +p) 


No 


[37] 


Scaling 


P 


- 


2/(3 +p) 


No 


[43] 


Scaling 


1 


1 


1/2 


No 




Shell 


1 


1 


0.5 


No 




Numerical 


1 


~ 1 


0.4 


No 


[3^ 


Scaling 


P 


1 


1/2 


No 




Scaling 


P 


1/2 


1/2 


Yes 


m 


Shell 


2 


- 


0.25 


No 


[36] 


Shell* 


P 


- 


1/(3 +p) 


No 


[38] 


EDQNM 


2 


~6/5 


~2/5 


No 


[39] 


Numerical 


2 


0.5 


- 


Yes 


[40] 


Numerical 


4 


0.7 


0.4 


No 




Numerical 


4 


1.1 


0.5 


Yes 


[53,5^ 


Numerical 


2 


1.05 


- 


No 




Numerical 


4 


0.5 


- 


Yes 


[47] 


Numerical 


4 


0.5 


0.5 


Yes 


[48] 


Mean-field 


P 


2(l+p)/(3 + p) 


2/(3 +p) 


No 




Mean-field 


P 


2/3 


2/3 


Yes 


[49] 


Numerical 


4 


~ 1 


~ 1/2 


No 


[50] 


Numerical 


4 


~ 1 


~ 1/2 


No 




Numerical 


4 


~2/3 


~2/3 


Yes 


[5^ 


Numerical 


4 


0.9 


~ 1/2 


No 




Numerical 


4 


-2/3 


-2/3 


Yes 



where pi is the matter density of ions and ain is the momentum transfer rate due to ion- neutral colhsions 53|, [54 1 . 

Although the form of the dissipation term is different before and after recombination, it is inessential to our 
discussion, since we work in the hypothesis of large kinetic Reynolds numbers, namely assuming that dissipation is 
negligible. [After recombination, the condition of tight- coupling is equivalent to assume that the dissipative term in 
Eq. (|138p is negligible in the Navier-Stokes equation [5J, [SJ] .] 

Proceeding as in section IVb, we can decompose the magnetic field according to Eq. (157)) . In this case, the quadratic 
term v • Vv can be approximated by v • Vv ~ Vq • Vv. This term in the Navier-Stokes equation can be neglected if the 
condition r|v|/|vo| ^ 1 is satisfied. In fact, comparing it with the Lorentz force, we have |v • Vv|/|Fi| — r|v|/|vo|. 
Since we are assuming |v| <C |vo|,J;he above condition is certainly valid if the F ratio does not assume very large values 
(and this is the case studied in j48| and in this paper). The approximation just discussed corresponds to assume that 
the small-scale components of the velocity field do not play any role in the development of MHD turbulence, at least 
on large scales. This means, in turn, that the Lorentz force acting on the charged particles of the fluid is responsible 
for the development of turbulence on large scales. In the limit of vanishing dissipation or high Reynolds numbers (see 
below), we have form Eq. ([31) that dtv ~ F^. 

To further simplify the Navier-Stokes equation, we write dt — I/t^, givirig v ~ TdF^, where r^ is the fluid-response 
time to the Lorentz force. This sort of "drag time' was first introduced in [43|, and its explicit expression as a function 
of time was derived in [48|, r^ = Td{0){l + r/700), where Td{0) = r(0) Toddy (0), and 

r(o) t 



— lim . 

t-^00 r(i) Teddy (i) 



(139) 



The quantity 700 is a constant since the results of numerical integration of turbulent MHD equations give, asymptoti 
cally, r{t) — 1 in the non-helical case and T{t) ~ constant (eventually different from unity) in the helical case 53|, |5J], 
while scaling arguments indicate that Toddy (0 ^ ^ 
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Consequently, the turbulent MHD equations are, formally, the same as the dragged MHD equations and then the 
master equation for the quantity ^diss , which rules then the evolution of the magnetic energy and correlation length is 

Looking at Eqs. (|140p and (|551) . we see that the substitutions S{0) — 1, 7(0) — >■ 700, and a = — 1 in the solutions found 
for the dragged case (see sections Vc, Vd, and Ve), give directly the solutions for the magnetic energy and correlation 
length in freely-decaying MHD turbulence. 

Vlllb. Evolution laws in radiation-dominated universe 

Non-helical case- Using the results in the above subsection and taking into account Eq. pi2p . we easily find the 
evolution laws of the magnetic field in radiation-dominated universe: 

RjBm - I ,,3 (iv,..)-^ (f )""" ,R>Ri. 

for the comoving magnetic field strength, and 

i£^^/'' ^ R^1^R^R^^ 

1/2 
for the comoving magnetic correlation length, where Vi = T^ B{Ri) [B{Ri) should be replaced by 

B[Ri)/ ^J p{Ri) + P{Ri) in the notation in which the quantity p + P is not taken to be 1, as already noted in section 
Vh]. In the above equations, 

^ c^ n^,{N,v,r\ (143) 

and K15, K16, K17 are given in Appendix A. In Eq. (|143p . we have assumed that Kij{NiVi)~^ > 1. In the opposite case, 
we have approximatively Ri ^ Ri. 

Helical case- In the helical case, we find 

1, Ri ^ R ^ Rl-: 



R'B{R) ^ j ^^^(M.,,A'mfjLr^ 



R2R,R,-i ■^..(^»^0-^(fp"^ Rl<R<R2, (144) 

R,B(R,) I _^ V ^^ _ xi/6 , X -1/3 



for the comoving magnetic field strength, and 

{1, ^ Rr<R<Ri, 

«i6 {N,v,)y('+p) (^)^ , Ri<R<R'2, (145) 

1/, / \-l/3 / x2/3 

Ki9 h]/\N,v,)^/^ (lnA^««.f ) (f ) ,R> R'2. 



for the comoving magnetic correlation length, where 



R2 H20 ,-2q/3 T- . -9„/-3.1 9/3 Jli, K91 . _9„/^ r. ,_9^/-!l29/3 

^^ fl T 



Ri NiV, 



's 



Inh-^^'^" 



-9/3 R!^ _ _K2^ 2g/3 



Ri NiVi 



In/j-^^/^ 



(146) 



with q — 3(3 +p)/(4p), and Kig, K19, K20, 1^21 are given in Appendix A. In the above equations, we have assumed that 
{NiViY'"^ hg "^ is a quantity much greater than unity. 
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VIIIc. Evolution laws in matter-dominated universe 

Non-helical case- In the case of a matter-dominated universe, we find 



for the comoving magnetic field strength, and 

UR) f^' R^<R<Ri, 



UR^) I ^23(^.«.)^/(^+^) (inf l'^'' ,R>Ri, ^^^^^ 



-R. 

for the comoving magnetic correlation length, where we have taken into account Eq. (I12ip . In the above equations, 

and K22, '«23, K24 are given in Appendix A. 

Helical case- In the helical case, and neglecting loglog factors, we find 

1 



1^22 ' 



T,^R(R^-\ ■^^AN^v,r^.(ln^J'''\ Rl<R<R2, (150) 

Ri^yR^) I ,,, ^ . ^ .-1/3 

'^25 /il/'(A^.«.)~^/' (in f j , i? > i?2, 

for the comoving magnetic field strength, and 

{1, ^ Rr<R<Ru 

^23 {N,v,f/(^+P^ (in f )^ , i?i < i? < i?^, (151) 

1/, / \2/3 

«:26/^s^'(W^«.)'/'(lnf) ,i?>i?^. 

for the comoving magnetic correlation length. Here, 



"-2 1 



i?,, \N,v,, ^ r R. \Nm, 



^M^h-,^'"]^ ^-exp ^/.-«/^), (152) 



and K25, '^26 J '^27, '«28 are given in Appendix A. The above equations are valid in the limit Hb ^ 1- 

Vllld. Exiting the turbulence phase 

The turbulence phase prior to recombination is characterized by high values of the kinetic Reynolds number. 
Re = vl/v 3> 1 where, as usual, v and I are the typical velocity and length scale of the fiuid motion. Since v, I, and v 
are function of time in an expanding universe, there could be a time when the system, initially in a turbulent phase, 
get a transition to a viscous phase characterized by Re < 1. 

After recombination, the kinetic Reynolds number, known in this case as "ambipolar" Reynolds number, has the 
form Rcamb = vlctinXe/ B^, whcrc Xe is the (constant) ionization fraction after recombination, and B is the typical 
magnetic field strength. Turbulence operates if Rcamb ^ 1, a condition that, although initially satisfied, can be (and 
indeed is) successively violated. 

Kinetic and ambipolar Reynolds numbers- In order to study the evolution of the kinetic and ambipolar Reynolds 
numbers, let as define them accurately as 

Re = , Rcamb = 55 ■ (153) 
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In a flat Friedmann universe, the kinematic viscosity re-scales as |34l.l53l.l54| 



\ R ^ly, radiation-dominated era, /,,,x 

iy^iy=< ,' (154) 

I R ' ly, matter-dominated era, 

where i/ as a function of the temperature is given by [38| i^(T) ~ [aomlog(l/Q^cm) T]^^ for T 3> m^ and i^iT) ~ 
{nf,aT)~^ for T ^ ?7ie, with acm being the fine structure constant, n^ oc T^ the electron density, and ax the Thompson 



cross section. The momentum transfer rate ain evolves with the expansion parameter as a-m oc R~^ 53l. |54||. so that 
in "tilde" variables we have ain = R^^^a^n oc R~^^^. 

Radiation- dominated universe.- In radiation-dominated universe, the asymptotic expansion of the kinetic Reynolds 
number is 

Re{R)o,h%, ^''>"^' (155) 

i? 3+p T <$: me. 



in the non-helical case, and 



^"^^^""l (lni?)V3i?-5/3^ r«m„ ^^^^^ 



in the helical case. 

Matter-dominated universe- In a matter-dominated universe, instead, we find 



i-p 



\ (lni?)^i?-5/2^ r<me. 



and 



^'^'^^"\(lni?)V3i?-5/2^ y<<^^^ (158) 

in the non-helical and helical cases, respectively. 

Regarding the ambipolar Reynolds number, we have Rcamb = Vi-msCvO-inXe/ B'^^^ in expanding universe. Using the 
fact that in a turbulent phase -Dims — tB^^^/S,b (see section Villa), we get Rcamb — Tai^Xe- Since f oc Ini? in a 
matter-dominated universe, we finally get 

Reamb(i?)oc(lni?)i?-3/2^ (159) 

in both non- helical and helical cases. 

Excepted the case of helical magnetic fields in radiation-dominated universe, we see that the kinetic and ambipolar 
Reynolds numbers are decreasing function of time (we remember that we are assuming p > \). Hence, there will exist 
two times icxit,i and toxit,2 when Re(ioxit,i) = 1 and Reamb(^cxit,2) = 1, thus defining the end of turbulence. '^ 

IX. ON THE CUT-OFF OF THE INITIAL POWER SPECTRUM 

,2 

In section Va, in order to get finite results for the magnetic energy, we introduced a Gaussian cut-off, oc e^ , for 
the initial magnetic energy spectrum. In the literature (see, for example, 39, 40, 47]), such a cut-off is sometimes 
associated to the initial configuration of the magnetic field itself, instead of its energy spectrum. This results in a e^ 
cut-off, sharper than the Gaussian one. Therefore, it is of some interest to consider the initial magnetic spectrum 

8B{k, 0) = Asfc^ exp(-2fc44), (160) 



^ Also for (maximally) helical fields there will be certainly a time when turbulence terminates. This because, after e+e~ annihilation 
around T^+^- ~ IMeV [59J], the kinematic viscosity increases by a factor of order 10^*^, namely of order of the photon-to-baryon ratio [59J. 
Accordingly, somewhere below T^j+j,-, the kinetic Reynolds number drops below unity. 
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instead of Eq. (|70] 



cut-off: Non-helical case- Following the same arguments in section Vb, we find for the magnetic energy and 



correlation length 



^-W-^-(«)^<4^'^'% 



CbW-CbCO) 



r( 4 ) ^\4'2^ 2 



r(^)t/fi±^,i% 



(161) 
(162) 



instead of Eq. (f77| and (f78| . respectively. In the limit of large Cdiss, which corresponds to large values of the time, we 
have the asymptotic expressions 



^,„ = ^,,o,^^rfi±£)ci'«'. 



1 r(^) 



(163) 
(164) 



where we used the asymptotic expansion of the confluent hypergeometric function U{a, b, z), namely U{a, h, z) ~ z~°' 
for z ^ 1 [56|, valid in the case a > 0. The master (differential) equation for ^diss is given by Eq. (|55)) . with the 
formal expression of (5(0) given by Eq. (|85l) . The solution of this master equation gives Eqs. (iQlT) and (1921) . with the 
constants C2 and C3 replaced by 



C2 



C3 



2(i+P)/-t f2,+p 



1 r(g±£) 

2l/4p^3+p\ 



4 

2(p-2)/2 



2V2r(3±£) 



C2, 



(i+p)r(i±i^V^' + ^ 



C3, 



respectively. Also, we find straightforwardly that the expression of 5(0) as a function of the index p is 

1/2 



m 



3V2r(i±H) 



r(^) 



1 



(165) 



(166) 



(167) 



instead of Eq. (|105l) . 



e cut-off: Helical case- Following again the analysis in section Vb, we find that in the helical case, the expressions 



for the magnetic energy, helicity, and correlation length, take the form 



Eb{0) 

HB{r) 
Hb{0) 

es(0) " 



Cl3(l 



/■p 

Set 



Ci4 



Cl5 



"^2p/2r[(i+p)/2]c]+f 

i + hB v^ cr' 

2 1^2 

Cb Cq 



exp 



exp 



Sq 
^Cdiss 

Sq 
2 Cdiss 



(168) 
(169) 
(170) 



in the limits Cdiss S> 1 and Ca /Cdiss 3> 1. Note that they differ from the case of a Gaussian cut-off for the introduction 
of the coefficients Ci3,ci4,ci5 (given in Appendix A), not present in Eqs. (I186p . (|187l) . and (|188l) . Therefore, following 
the analysis in the Appendix B, we find that the expression for the magnetic energy and correlation length are given, 
respectively, by Eqs. (|97| and (|98|) . with the replacements 



C7 



C8 



r(i±H)r(|) 




r(i±H)r(H) 



-2/3 



(i + p)r(i±£) 
V2r(^) 

(i+p)r(i±£) 



-1/6 



C7, 



-,1/6 



C8- 



(171) 
(172) 
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Other cut-offs- We conclude this section by observing that very similar results are obtained if the cut-ofF of the initial 
spectrum is even sharper than the e~ one, as for example in the case 

£B{k,0)^h^^'^ ^<^' (173) 




In general, what we observe is that changing the form of the cut-off, just changes the form of the coefficients C2, C3, 
C7, Cg, and the expression for S{0) as a function of p. This, in turn, changes the form of the k's coefficients which 
enter in the expression for the magnetic field intensity and correlation length in the case of interest of an expanding 
universe. However, the power-law behaviors, and logarithmic corrections in the helical case, do not depend on the 
form of the cut-offi 

Finally, we note that the authors in [50], as already explained in section Va, found that soon after a first-order 
phase transition, the magnetic spectrum has the form of a Batchelor spectrum £b oc fc* at small wavenumbers. Also, 
they found that during a turbulent phase, a Kolmogorov spectrum of the type £b oc k~^'^ develops at intermediate 
wavenumbers, kB < k < fcdiss, where fc^ = Stt/^b and fcdiss = 27r/-y/77t, with ^/rjt being the dissipation scale. This 
means that, if a dragged phase follows a turbulent phase (as it could be in a cosmological context 53|, l5J|), a more 
realistic form of the initial magnetic spectrum in dragged phase could be given by 

{Aifc*, k < kB (Batchelor spectrum), 

X2k~^^^, kB < k < kdiss (Kolmogorov spectrum) , (174) 

A3exp(— 2fc^//c^jgg), fc > fcdiss (dissipative cut-off), 

where Ai, A2, and A3 are constants, fc^ = 27r/^B(ii), ^diss = ^t^ I \fWi^ and ti the initial time. So, in this case the 
cut-off is represented by a Kolmogorov spectrum up to the initial dissipative scale, 27r/fcdissj followed by a Gaussian 
cut-off due to resistivity. 

In the light of the above discussions, however, we expect that even in this more realistic and complicated case, 
the main results we found, namely the power laws for the magnetic intensity and correlation length (corrected by 
logarithmic factors in the helical case), do not change. 

X. DISCUSSION 

In this section, we discuss our results and write down the final and simplified expressions for the magnetic field in- 
tensity and correlation length in a Friedmann universe, which are relevant when studying the evolution of a primordial, 
phase-transition-generated, cosmic magnetic field. 

Simplified evolution laws— We saw, in the previous section, that the n coefficients entering in the expressions for 
the magnetic field and correlation length depend on the choice of the cut-off of the initial magnetic spectrum. In the 
case of magnetic fields generated in primeval phase transitions, the exact form of this cut-off is not precisely known, 
and this introduces a factor of arbitrariness in the evolution laws of B and ^b ■ Fortunately, since the k coefficients are 
all of order unity (see section IX and Appendix A) , this factor of arbitrariness can be safely neglected when studying 
the evolution of a phase-transition-generated magnetic field. 

If on the one hand the k coefficients depend on the initial cut-off, on the other hand the logarithmic factors we 
found in the case of helical magnetic fields are a general consequence of the fact that the magnetic helicity is a 
quasi-conserved quantity in magnetohydrodynamics (see Appendix B). However, these factors introduce just a small 
correction in the power-law evolution of B and ^b and then, in first approximation, can be neglected. 

According to the above discussion and given an initial magnetic field with spectrum of the form £B{k, ti) ex fc^, the 
(physical) magnetic field intensity, B, and the comoving magnetic correlation length, ^s, evolve as a function of the 
temperature T approximatively as 



B 
B, 



f(ft) 



Ti<T<Ti, (quiescent phase), 
Y-j , Ti<T< T2, (selective decay phase), (175) 

h]^^{N,v^y^'■°'^ (^X'' , T>T2, (inverse cascade phase). 
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TABLE II: The exponents Qi{p), Q2{p), Qsip), and Q4,{p) in the evolution laws p75|) and (|176[) of a primordial magnetic field in 
dragged and turbulent magnetohydrodynamic phases, in radiation (RD) and matter (MD) eras. For the dragged phase, it is 
indicated the free-streaming species which determines the drag coefficient. The parameter p is the index of the initial magnetic 
power-law spectrum, £B(k,ti) oc k^ . 



Phase 


Era 


Streaming particle 


qAp) 


f>2(p) 


Q-iij>) 


Qiij>) 


Dragged MHD 


RD 

MD 

RD 

MD 


neutrino 
photon 
photon 

hydrogen 


1+P 

2(3+p) 
1+P 

2(3+p) 
1+P 

2(3+p) 
1+P 

2(3+p) 
1+P 
3+p 
1+P 
3+p 


17+9p 
2(3+p) 
4(2+p) 

3+p 
33+17p 
4(3+p) 
39+23p 
4(3+p) 
7+3p 
3+p 

2 


1 

3+p 

1 
3+p 

1 
3+p 

1 
3+p 

2 
3+p 

2 
3+p 


5 
3+p 

4 
3+p 

9 




2(3+p) 
15 


Turbulent MHD 


2(3+p) 

2 

3+p 





and 






1, 



(*) 



Biiv) 



Tt<T <Ti, (quiescent phase), 

^1 ^T <T2, (selective decay phase). 



(176) 



h]i'{N,v,r^(^^ (^) 



23(0) 



,T>T2 



(inverse cascade phase). 



respectively, where Bi — B{Ti) and ^B,i — £,B{Ti), and we took R ^ T^^ [neglecting small corrections due to (7*^s(T)]. 
We remember that Ni is the initial number of magnetic domains per horizon, namely Ni = dH,i/£,B,phys,i, where dH,i 
and ^B,phys.i are the length of the Hubble horizon and the physical magnetic correlation length at the initial time. 
The value of the bulk velocity at the onset of the particular regime, Vi, is not an independent parameter in both 
dragged phase and turbulent non- helical phase, being related to Bi, ^B.phys,i, and at by the relations 



ai5B,phya,i 

Bi, 



(dragged phase), 

(turbulent phase - non-helical case) 



(177) 



In equations (|175l) and (|176p . we quoted the expressions for B and ^^ in the helical case. In the non- helical case, 
those expressions are still valid, the inverse cascade phase being absent. The transition temperatures Ti and T2 can 
be easily found by matching the expressions of B (or, which is the same, those for S.b) in two consecutive phases. 

We summarize our results on the evolution laws for the magnetic field and correlation length in Table 2. 

Comparison with the Banerjee and Jedamzik 's results- Let us now compare the above results with those of [53l |54| , 
which are, to our knowledge, the only ones in the literature which consider other than the case of freely-decaying 
magnetic fields, also the case of magnetic fields in dragged phase. To this end, we introduce the so-called Alfven eddy 
turnover time, r^, as ta — CB,phys/B. We can then write the quantity NiVi as 



N,v,, 



{ti/TA.i){Td,i/TA,i), 

ry\u/rA,), 



dragged MHD, 
turbulent MHD, 



(178) 



where we remember that in turbulent regime, Tj is of order unity for the non-helical case, but can be different from 
one in the helical case. Banerjee and Jedamzik 53|,|54| assume (without a full justification) that the initial correlation 
length is not a free parameter but it is determined, at the initial (cosmic) time ti, by the equality of the Hubble rate 
H ~ t~^ and the Alfven eddy turnover rate 1/ta- Namely, in their model we have 



ti ^ TA,: 



In this case, and taking F^ of order of unity (as in [53|, l5^), we get 



N,,v,, 




dragged MHD, 
turbulent MHD. 



(179) 



(180) 



With the choice in Eq. (I179|) . it is easy to check that our formulae agree with those of [53|, |5J| [with the exception of 
Eqs. (62) and (63) of t53..54|. in which the exponents 3/(2 -I- n) and 3n/(2 -I- n) should be 4/(2 -I- n) and 4n/(2 -I- n), 
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respectively.] For example, for non-helical dragged MHD in the case of neutrino free streaming we have, excluding 
factors of order unity, 

l+p 17+9p 

/ 1 \"2(3+p) /T^X 2(3+p) 

^-(^) - ^-(^^) {cf^T (I) "" ' "^ >^ "^^ ^^''^ 

where Gp is the Fermi constant, mpi is the Planck mass, and we used the facts that ti ^ mpi/Tf (in radiation- 
dominated e ra) and Ui ~ GpTf (for the neutrino case). Equations ()18ip and (|182p agree, respectively, with Eqs. (48) 
and (47) of ^. [Let us observe that in the notation of Q, p = n - 1 and B{T) ex r{TY''^T'^.] 

Inertial range- According to the "Kolmogorov hypothesis" (see, e.g., ^23||), the cascade of energy in fc-space is a 
quasi-local process independent on the particular scale considered (although this is strictly possible only for scales 
larger than the dissipation scale and smaller than the integral scale). The fc-space interval where the Kolmogorov 
hypothesis may apply is called "inertial range". Here, the magnetic spectrum is expected to decays as £B{k,t) oc fc^, 
with /3 being negative. The value of /3 could in principle depend on the particular phase (dragged or turbulent) and be 
different for the non-helical and helical cases. Indeed, the pure Kolmogorov spectrum, (3 — —5/3, has been observed 



in [50| for the turbulent non-helical case, while [53|,|5J1 found /3 ~ —2 both in helical and non-helical dragged case. 

In our case, we do not find any inertial range where the magnetic spectrum decays as a power law. This is probably 
to ascribe to mean-field approximation, which neglect the small-scale fluctuating part of the velocity and magnetic 
fields (see section IVb). Is the (quasi-linear) interaction of these small-scale fields with the corresponding large-scale 
ones which "opens" this particular range in fc-space. However, if the inertial range is sufficiently narrow, we do not 



expect significant modifications to the evolution laws we found for B e ^b- Indeed and interesting enough, Ref. [53|, [54 1 
found that the inertial range is almost absent in the non-helical turbulent case, and that on intermediate scales the 
magnetic spectrum is more consistent with an exponential law than a power law. This seems to be in agreement with 
the results we found in mean-field approximation, namely with Eq. (j94]). 



XI. CONCLUSIONS AND OUTLOOK 

Conclusions- The presence of large-scale magnetic fields in the universe is still an open and unsolved mystery of 
modern cosmology. A plethora of mechanisms able to generate cosmic magnetic fields in the early universe, as for 
example during inflation or primeval phase transitions, or during more recent eras by astrophysical mechanisms have 
been proposed since the first attempt put forward by Harrison in 1970 [ll|. 

If on the one hand, astrophysical mechanisms seem to be ruled out by recent observations of large-scale magnetic 



fields in cosmic voids 17|, on the other hand inflationary mechanisms are able, at least in principle, to explain 
the magnetization of the universe. However, these latter mechanisms repose on the use of nonstandard physics, an 
exception being the mechanism recently proposed in [22|. Nevertheless, if the scale of inflation is considerably below 
lO^^GeV, also this mechanism fails to explain galactic and galaxy cluster magnetic fields. 

If, instead, presently-observed magnetic fields originate from cosmological phase transitions (such as QCD or elec- 
troweak phase transitions), another unsolved question arise: how have these magnetic fields evolved from the time of 
they generation until today? The answer to this question can be given only in the framework of magnetohydrody- 
namics. 

In this paper, we have studied the evolution of phase-transition-generated magnetic fields coupled the to primeval 
plasma by solving, analytically, the magnetohydrodynamic equations in mean-field approximation. The reduction 
of full MHD equations to simpler equations due to mean-field approximation causes a lost of information about 
the transfer of magnetic energy at intermediate scales (corresponding to the inertial range). Nevertheless, the main 
characteristics of an evolving magnetic field are preserved. 

In particular, we have analyzed the decay of primordial magnetic fields, both in radiation and matter eras, in the 
two regimes which are relevant in a cosmological context, namely the turbulent and viscous free-streaming regimes. 
During a viscous free-streaming phase, dissipation is ruled by a drag term in the Navier-Stokes equation, the drag 
term depending on the free-steaming particle species, namely neutrinos, photons or hydrogen atoms. 

Our results can be summarized as follows. 

If the initial magnetic field is not helical, then the magnetic helicity remains null during the evolution of the 
system. There is an initial phase in which the system is quiescent: magnetohydrodynamic effects do not operate 
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and the characteristic quantities of the system, such as energy and correlation of the magnetic field, remain almost 
constant. This phase persists for a period which is proportional to the initial eddy turnover time. After that, the 
evolution of the magnetic field proceeds through selective decay of magnetic modes, that is there is no direct transfer 
of magnetic energy from small scales to large scales, but simply modes with larger wavenumbers are dissipated faster 
than those whose wavenumbers are small. During this process, the magnetic correlation length grows, while the 
magnetic energy decays in time. The evolution laws for the magnetic field depend on the initial magnetic power 
spectrum and on the particular regime. They are summarized in Eqs. (|175[) - (|177p and Table II. 

In the helical case, the magnetic helicity is an almost conserved quantity. The system undergoes three different 
phases: a quiescent phase and a selective-decay phase in which the system evolves irrespective of the presence of mag- 
netic helicity, and an inverse-cascade phase. As in the case non-helical case, the quiescent phase ends approximatively 
after a period of time equal to the initial eddy turnover time. Then, the system enters into a selective-decay phase 
characterized by a magnetic field evolution similar to the non-helical one. This last phase ends due to the conservation 
of magnetic helicity, favoring an inverse cascade of the magnetic field during which the magnetic energy stored on 
small scales is partially transferred to larger scales. This process of energy redistribution weakly depends on the prop- 
erties of the initial magnetic field, so that the evolution laws of magnetic energy and correlation length do not depend 
on the index of the initial magnetic spectrum power law, although they are different in turbulent and free-streaming 
regimes [see Eqs. (|175p - (ll77p and Table II]. Moreover, the time when the system enters into the inverse-cascade regime 
depends both on the form of the initial magnetic spectrum and on the amount of initial magnetic helicity. 

Outlook- Our analytical evolution laws in the different regimes and for helical and non-helical magnetic fields. 



are in substantial agreement with the numerical results obtained in [53l. l54j. Nevertheless, it is important to stress 
that (i) there is no yet convergence in the literature on the results of numerical simulations of MHD turbulence 
in the case of non-helical fields, and that (ii) the only available numerical results for the dragged phase are those 



53|, |5Jj. Therefore, it is desirable, in order to have a firm understanding of the evolution of primordial magnetic 



fields, to solve this disagreement for turbulent non-helical fields and, at the same time, to have an independent (and. 



possibly, full numerical) confirmation of the Banerjee and Jedamzik results 53|, |5J| for the evolution of magnetic fields 
in dragged phase. Finally, (iii) we stress that the study in ^, |5J| on the properties of phase-transition-generated 
magnetic fields throughout the evolution of the universe assumed a linear dependence between the initial magnetic 



correlation length and the initial magnetic field intensity. It would be important to see if the conclusions of 53|, l54j , 
that phase-transition-generated magnetic fields may directly account for galactic and galaxy clusters magnetic fields, 
are modified or even invalidate by relaxing the above ansatz on the initial magnetic correlation length [69|. 

APPENDIX A 

In this appendix, we write down the expressions for the coefficients c; and Ki. They are, respectively, 

(3 + p)cI^(o)7(o) „ _ JIs + pKIMIiM!!'*' . _.-(i+p)/2 , _ 1/2 

C2 — C-^ 1 C.^ — ^1 ■! 

3 ^;:x^^^^„,^nMa^„^l/3 .^3 ^,1/3 



'° = 6(1^^^ '^^^1 6(1^^^ ^ '"^ = "1 '"^ = "1 



c. = 24[W0)]-^c„ C5-{fV! '^"!^°'' } ' ^a = {f P^^(0)[7(0)r/..(l-a)} 



C7 



I 1 ^^^ 

3J(0H7(0r ^^^^'"7 '^« = '^^^"' 



-Q 



2(l + p) 2(3 + 2p) 



Cl2 = ICb 



^ PS 2q/3 

^^ x/2(21n2+p) f (3+p)J(0)[7(0)]M — f 
Vln2+p I 6(1 -a) 



2(ll+3p)/4p/-l±p-\ 23(4+P)/4 
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and 



KO 



K4 



KS 



1+p 

12 



1/2 



1/6 



Cs(3+P) 



i+p 

'2(3+p) 



^, _ ^-2/(l+P) ^, _ 2(3+p)/(l-a)(l+p) 



6a/3 

Cb(3+p)3" 



(l-a)v/3(l+p) 



1+p 

2I3+PT 4(3+p) 

-2/(l+p) - 3a(l + p) 

, Kg = Kg " ', Kio = Kg ,^11 



6%/3 



1/6 



(_a)l/3 3-a/6^ 



K12 = Kji , Kl3 = (kh/ks) '', Kl4 = (K12/K9) '', Kl5 



-2q 



(3+p)c: 



1270 



'^16 — 1^15 



2/(l+p) 



^p)k15, K23 = 2-2/(3+p)ki 
K27 = (k25/k22)"^'', K28 = (k26/'«23)"^''- 



«22 = 2(i+p)/(=^+p)ki5, «:23 = 2-2/(3+P)«i6, K24 = ng^^^ ^ ^'"'''K K25 = 2'^'n^s, K26 = 2-2/3«:i9, 



(184) 



APPENDIX B 



In this appendix, we find the asymptotic solution of Eqs. (l551) - ([5^ . where Eb and Hb as a function of Cdiss and Ca 
are given in Eqs. ([77| and (|78|) . respectively. 

In the limit Ca/\/l + Cdiss ^ 1' ^^'^ magnetic energy and correlation length, Eqs. (f77|) and ([8T|). are given by 
Eq. ([Ml) and ^, respectively Inserting Eq. dHU) in Eq. ([83|), we get 



f]/-2 a2 a2 _i+p 



(185) 



Consequently, in 



whose solution, for very large initial magnetic Reynolds number, ReB(O) 3> 1, is given by Eq. 
the case Ca/\/l + Cjiss ^ 1' ^'^^ magnetic energy and correlation length are given by Eqs. (|9Tt and ((92l) . respectively, 
to wit, the same as in the non- helical case (see section IV a). In the limits Cdiss ^ 1 and Ca /Cdiss ^ 1, the expressions 
for the magnetic energy, helicity, and correlation length, Eqs. ([77]) . (|78p . and (|8T|) . take the form 



(l + /iB) 



2p/2r[(i+p)/2]ci+^^ 



C; 



p— 1 



i^B(r) 
£b(0) 

Hb{t) 1 + hB 

Hb{0) ^ hB 2(p-i)/2r(p/2) q:+'^ 

Cb(0) Cb Co ' 



^" exp ^" 



2Cd^, 



diss 



exp 



K 



2 

■diss 



(186) 
(187) 
(188) 



respectively, where we used the asymptotic expansion of the Kumnier confluent hypergeometric function ii^i(a, b; z) 
for z » 1 ^, iFi(a, 6; z) -- [r(fe)/r(a)] z'^^^'e^, valid in the case a > 0. From Eqs. p^ - ([TMl) . one directly obtains 
Eq. (|101[) . Because the magnetic helicity is quasi-conserved, Hb{t) ~ Hb{0), from Eq. (|187l) we get 



exp 



^ Cdiss 



2(p-i)/2r(p/2)Cd 



-l+2p 
iss 



l + hB 
Inserting the above equation in Eq. (|186p . we get 






Cb Co 



Eb{0) - 2 C^,,, 



(189) 



(190) 



For r ^ 1 and from Eq. (jl89p . we have Cq, ~ -\/2Cdiss[(2p— 1) In Cdiss + (1 ^p) hiCa]^^^- From the above expression we 
get, neglecting terms of order In (In Cdiss), 



Ca^(2plnCdiss)'/'C 



(191) 
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Inserting Eq. (IT9T1) in Eqs. p90)) and (IT88)) . we obtain 

Cs(0)-VpCb (InCdiss)!/^' ^ ^ 

respectively. In the limit of very large initial magnetic Reynolds number, and taking into account Eq. (|192p . equa- 
tion (1551) reads 



^ = ^ .[^SiOMOrhB r-^ Miif^ . (194) 

«T 6 \ Z (.diss 

Note that the above equation is valid only for t 3> 1. Introducing, for convenience, the reference time t^ such that 
T* :^ 1, the solution of Eq. (|194p is given by 

erfi [v/31nCdiss(r) ] - erfi [^Sln^Ur.) ] ^ \[^ ^ M ^^^^^^^^^^^^"^ (r^^ - ri^) , (195) 

V7r6V2 1 — a 

where erfi(z) is the imaginary error function 56]. By using the asymptotic expansion of the imaginary error function 
for z >• 1 ^56], erfi(z) ~ [l/(\/7rz)] e^ , we can write Eq. (|195p . for t ::^ t^,, as 

C„..(.)m/| ^(»)W°)l°"« f(i„a,..)./.,i.-.,/3^ ,130) 

[^ Z y 2 I — a J 

Neglecting terms of order In(lnT), we can write the above expression as 

[ z \ b a/I - a J 

which is indeed Eq. (|95|) . Now, inserting the above equation in Eqs. (J19ip . (jl92p . and (jl93p . and neglecting terms of 
order In(lnT), we get 

<.M = {f p'i(0)(7(0ri.fl(l-o)| (hirf'T"-"", (198) 

respectively, which are indeed Eqs. (|96|) . (|97| . and 



APPENDIX C 

In this appendix, we derive Eq. (|103p . namely the expression for the time tg when the inverse cascade begins. 
In section Vd, we defined it as the time when the maximum of the magnetic spectrum meets the initial spectrum, 
namely SBikmaxjtg) = fs(fc,nax,0), where fcmax is determined by the condition d£Bikma.x,t)/dk — 0. Introducing the 
quantity x = kf.B, the above two conditions read 

cosh(2a;,„axCQ) + ^b sinh(2a;,„axCQ) = exp(2a;f„axCdiss)' (201) 

2xmaxCa[/is cosh(2xmaxCa) + sinh(2xmaxCa)] + [v - 4a;^ax(l + Cdiss)] exp(2x^^^Ciss) = 0, (202) 

respectively, where Xmax = fcmax^s, and all quantities depending on the time are evaluated in t ~ Tg. In the limit 
Hb < 1 and 2x„iaxCa > 1, Eqs. ^UB) and ^U^ become 2x„iaxCa = ln2 + 22;^^^Ciss and 2x^ax(l + Ciss) = ln2 +p, 
respectively. In the limit CdissiTg) S> 1, these equations read 



2XinaxCa-21n2+p, (203) 

inaxSdiss 



2a^LxCdiss = ln2+p, (204) 



35 

respectively. From Eq. (I203p . we get that the condition 2a;,„axCa ^ 1 is indeed satisfied for p > 1. Moreover, since 
for T = T£ we have from Eqs. (|203l) - ((204)) that 

Cc 21n2+p , ^ 

^ < 1, (205) 



Cdiss ^2(ln2+p) 

we can use Eq. (j90| . valid for Ca/Cdiss ^ 1 and t ;» 1 (sec section Vd), evaluated at r = rg, to wit Cdiss = 
cyV|^""^/^^+^\ This equation is then valid for rg > 1. Using Eq. ^ in the hmits Hsit) = iJs(O), Ca/Cdiss < 1, 
and Cdiss ^ 1, we get Ca/Cdiss = /ibC-B c\' t|^ ^'^)/( ^'^2. Comparing this equation with Eq. (|205p . we find 

_ /^-^^ V2(21n2 +p) f (3 +p)J(0)[7(0)]'- ^'^f ^^' ^2^/3 
"^^-f^ Vln2+p 1 6(1^^^ / / ^^ ' ^^"^^ 

which is indeed Eq. (|103p . Since /is ^ 1, the conditions rg ^ 1 and Cdiss(''£) ^ 1 are satisfied. 

Finally and for the sake of completeness, we give the expression for /cmax- Taking into account Eq. ()204p . we have 
kmnxis = a^max = (In 2 + p) ^^^ [2ci (C12) p« ] ~^^^ ^B ■ Notc that k^nx^B is a small quantity in the hmit of small hs- 
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